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In order to investigate the systematics of the loop expansion in high temperature gauge 
theories beyond the leading order hard thermal loop (HTL) approximation, we calculate 
the two- loop electron proper self-energy E in high temperature QED. The two-loop bubble 
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diagram of E contains a linear infrared divergence. Even if regulated with a non-zero photon 



mass M of order of the Debye mass, this infrared sensitivity implies that the two-loop self- 
energy contributes terms to the fermion dispersion relation that are comparable to or even 
larger than the ncxt-to-lcading-order (NLO) contributions of the one-loop E. Additional 
evidence for the necessity of a systematic restructuring of the loop expansion comes from 
the explicit gauge parameter dependence of the fermion damping rate at both one and two- 
loops. The leading terms in the high temperature expansion of the two-loop self-energy for 



all topologies arise from an explicit hard-soft factorization pattern, in which one of the loop 
integrals is hard (p ~ T) , nested inside a second loop integral which is soft (0 < p < T for real 



parts; p ~ eT for imaginary parts). There are no hard-hard contributions to the two-loop 
E at leading order at high T. Provided the same factorization pattern holds for arbitrary I 
loops, the NLO high temperature contributions to the electron self-energy come from I — 1 
hard loops factorized with one soft loop integral. This hard-soft pattern is a necessary 
condition for the resummation over I to coincide with the one-loop self-energy calculated 
with HTL dressed propagators and vertices, and to yield the complete NLO correction to 
E at scales <~ eT, which is both infrared finite and gauge invariant. We employ spectral 
representations and the Gaudin method for evaluating finite temperature Matsubara sums, 
which facilitates the analysis of multi-loop diagrams at high T. 
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I. INTRODUCTION 

The systematics of the usual loop expansion in powers of the weak coupling e 2 = 1/137 <C 1 is well 
understood in quantum electrodynamics (QED), in perturbation theory about the vacuum state. However, 
when one considers extreme states of matter at high temperatures or densities, straightforward perturbation 
theory is of little use, even in a weakly coupled gauge theory like QED. At very high temperatures (T 3> m e ) 
one signal of the failure of the usual loop expansion is the appearance of infrared singularities of the form 
(eT/ta) 2i ~ 2 , at I loop order, where w is a characteristic frequency, such as that of an external line in 
a self-energy diagram. The appearance of such infrared singularities implies that higher order terms in 
the perturbative series in e 2 are unsuppressed for low energies or momentum scales, and hence the usual 
perturbative loop expansion fails for oj < eT. 

The singular infrared behavior of perturbation theory at finite temperature can be seen already at 
the one-loop level. The leading order one-loop singularities arise from the ultraviolet or hard region of 
the momentum integral, whose power law growth is cut off only by the Bose-Einstein or Fermi-Dirac 
thermal distribution, so that the loop momentum becomes of order of the temperature T. These leading 
contributions in the high temperature expansion of one-loop self-energies and vertices are the hard thermal 
loops (HTL's) [1], studied extensively in the 1990's in both QED and QCD (see [2-5] for reviews). It was 
found that the leading order (LO) HTL self-energies and proper vertices enjoy some remarkable properties 
[6-8] , among which is that they obey the same Ward identities as those of the exact theory. This enabled 
Braaten and Pisarski to propose an HTL resummation program [6, 8], in which the LO HTL propagators 
and vertices are used instead of the bare ones in any diagram in which the loop integration momentum 
becomes soft (~ eT <C T). This resummation algorithm was successfully applied to a number of problems, 
chief among them the first determination of the damping rate of gluons at rest in the high temperature 
QCD plasma, which is finite and free of any gauge ambiguities [9]. 

The Ward identities obeyed by the lowest order HTL self-energies and vertices, and the success of HTL 
resummation in producing gauge invariant physical quantities suggests it should be possible to restructure 
the perturbative loop expansion of gauge theories at high temperatures or densities in a systematic way. 
To do so, however, requires a more complete understanding of the infrared singularities of the subleading 
terms in the high temperature expansion at higher loops. The subleading or next-to-leading order (NLO) 
terms in the high temperature expansion of the one-loop level fermion self-energy have been computed 
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[10, 11], but there have been only a few partial calculations of various quantities in different kinematic 
limits at two-loop order [12-16]. As a consequence, both the nature of the subleading infrared singularities 
of perturbation theory which HTL resummation includes and those which it does not include, and the 
general systematics of the reordering of the loop expansion at high temperatures beyond HTL's, if it 
exists, has remained relatively unexplored. 

The promising beginning of HTL resummation of infrared singularities, and the possibility of extending 
its success to higher orders in a systematic way is especially interesting in gauge theories. More than just 
a formal question about the nature of perturbation theory at high temperatures, the lack of systematic 
expansion methods that can capture both collective behavior and infrared physics in hot or dense matter, 
and also preserve gauge invariance and the UV renormalization structure of the underlying microscopic 
theory is a serious obstacle to physical applications in processes from the early universe to heavy-ion 
collisions, both in and out of equilibrium. 

Motivated by these considerations, we present in this paper a calculation of the two-loop electron 
self-energy £ in high temperature QED. Our approach is simply to calculate the two-loop self-energy 
in the high temperature limit, systematically extracting the leading terms at high T, without making 
any assumptions about the relative importance or unimportance of different contributions or momentum 
ranges, and without assuming a priori how they should be resummed. Once these leading terms in the high 
temperature expansion of the two-loop self-energy are known, they may be compared then with the known 
NLO terms of the one-loop self-energy [10, 11]. By explicit calculation we shall find that the leading terms 
in the high temperature expansion of the two- loop self-energy are of the same size, or even larger than the 
NLO terms of the one-loop Si when evaluated on the fermion mass shell. This is a direct demonstration 
of the need to take two and higher loops into account in order to obtain the full NLO result. 

Dissecting the explicit two-loop self-energy calculation also allows for a more complete understanding 
of the relative importance of different momentum ranges in the loop integrations. It is already well known 
that the LO HTL result at one-loop order results from the strictly hard region of the loop integration 
where p ~ T. The NLO terms of the one-loop £ arise from the region of the loop momentum integral 
which is parametrically less than T. For the case of the imaginary part of S, the relevant momentum 
scale is p ~ eT <C T, while the real part of the self-energy receives contributions form the entire range of 
loop momentum < p < T, and give rise to ln(T/w) weak dependence upon the temperature. Although 
momenta up to T do play a role in the real parts, unlike the purely hard contributions which grow as a power 
of T, logarithms are a signal of sensitivity to the low momentum, infrared region of the loop integration as 
well. Following the convention in the field [1, 6], we term all such contributions where internal momenta of 
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order eT or less play an important role, soft loop integrals. We shall see that the leading order terms at two 
loops which are comparable to the NLO terms at one loop all come from momentum integrations where 
one of the loops is hard, while the other is soft, in the sense just described. Thus these comparable terms 
coming from one- and two-loop diagrams have in common that they arise from exactly one soft internal 
momentum loop. 

An additional signal of the need for a reordering of perturbation theory comes from the gauge dependence 
of the one- loop self-energy, even when evaluated on the fermion mass shell. In the LO HTL limit, i.e. when 
the loop momentum is hard, the one-loop electron dispersion relation is independent of the gauge fixing 
parameter £. However at NLO in the high temperature expansion coming from the soft region of the 
loop momentum integral, the one-loop dispersion relation becomes dependent upon £ [10, 11]. This is 
sufficient proof that the NLO one loop calculation is incomplete and physically meaningless by itself, and 
a systematic resummation method to obtain gauge invariant results for physical quantities such as the 
fermion damping rate at order e 2 T is necessary. Our calculation will show that this gauge dependence is 
not cured at two-loops, and a full resummation of higher loops is necessary in order to arrive at the full 
NLO gauge invariant decay rate of an on-shell fermion in the plasma of order e 2 T. 

When higher numbers of loops are considered, we encounter a new phenomenon not present at one 
loop. At two-loop order in addition to gauge parameter dependence of the on-shell dispersion relation, 
we encounter uncanceled true infrared and/or collinear divergences, not regulated even at finite values 
of the soft external energy u, or momentum q. In particular, a linear infrared divergence appears in 
the bubble diagram Fig. 2 of the two-loop fermion self-energy, arising from the double photon massless 
propagator pole in this diagram. Quite apart from dependence upon the gauge parameter £, the appearance 
of such unregulated divergences from massless, unscreened photons is itself also sufficient evidence for the 
breakdown of the loop expansion at high temperatures, and the need for a resummation of diagrams with 
still larger numbers of loops. 

In order to regulate the divergences in the two-loop fermion self-energy, we introduce a photon mass M 
in the evaluation of the bubble diagram of Fig. 2. In the usual perturbative expansion, the formal limit 
M — > should be taken. Since this limit does not exist at two-loop order, we may anticipate the result that 
after HTL resummation of the photon self-energy the photon acquires a non-zero Debye screening mass 
of order eT. Taking M ~ eT allows us to characterize the order of all the two-loop contributions to S, 
and compare finite contributions from different loop orders in a meaningful way. In other words, when the 
formal photon regulator mass is replaced by the physical Debye screening mass, by a partial resummation 
of the photon propagator for soft momenta, there are neither infrared nor collinear divergences in the 
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two-loop dressed self-energy, and all NLO contributions to the fermion self-energy are estimated correctly 
in powers of e. The photon mass M thus plays a dual role in our considerations, first as a regulator of 
infrared divergences in bare loop diagrams, and second, as a parameter whose finite value of order eT after 
resummation of higher loop diagrams enables us to estimate the finite magnitude of resummed amplitudes. 

Both the gauge parameter dependence and the explicit infrared divergences as M — > in the bare loop 
expansion can be cured only by systematically including the contributions from higher numbers of loops 
£ > 2, despite the fact that these contain additional explicit powers of the coupling. This HTL improved 
resummation leads finally to the complete order e 2 T corrections to the fermion dispersion relation in the 
plasma for both its real and imaginary parts, compared to the LO result for the fermion thermal mass, 
eT/y/8 which is purely real. Thus the reordering of perturbation theory in the plasma is decidedly non- 
analytic in a = e 2 /Ait. 

The structure of the explicit two-loop calculation of £ shows that the leading contributions at high 
temperatures arise entirely from momentum integrals which factorize into one momentum integral which is 
hard (p ~ T) , and a second one which is soft (in the sense defined above) . Overlapping two- loop momentum 
integrals which violate this hard-soft pattern can be shown explicitly to be suppressed by additional powers 
of e, thus justifying their neglect, in 0{e 2 T) calculation of the self-energy and damping rate. In particular, 
the detailed analysis of the two-loop self-energy shows that there are no contributions at this order coming 
from both loops hard, which naive power counting would apparently allow [3, 8]. That such hard-hard 
contributions from two-loop diagrams at this order in e 2 T are completely absent is in accordance with the 
an earlier result by Schulz for the two-loop gluon polarization in QCD at next to leading order [14]. A 
fact apparently not widely recognized is that the contribution at two-loop order of only hard-soft integrals 
to terms of order e 2 T in the fermion self-energy and the complete absence of hard-hard contributions 
is exactly what is required by the HTL resummation program, when HTL propagators and vertices are 
re-expanded in primitive bare loops of perturbation theory, a point which we elaborate in some detail in 
Sec. V. Thus one result of our detailed evaluation of the two-loop self-energy is a refinement of naive 
power counting rules, and a clarification of the hard-soft factorization structure of the high temperature 
expansion of perturbation theory, which underpins resummation of one-loop diagrams with HTL dressed 
propagators and vertices. 

The factorization pattern should be expected to continue at £ > 2 loops as well, namely all the order 
e 2 T corrections to the fermion dispersion relation should arise from £ — 1 hard loops inside only one 
soft loop. This HTL reordering of the ordinary loop expansion is necessary to collect all the infrared 
sensitive terms of relative order (eT/u) 2t ~ 2 from £ loops, which make equally important contributions to 



6 



the fermion dispersion relation and damping rate, when uj ~ eT. When all these subleading contributions 
are summed over all t loops by the HTL resummation of propagators and vertices in an HTL dressed 
one-loop diagram, the gauge parameter dependence at subleading e 2 T order is eliminated. The curing of 
both the infrared divergences, and the gauge parameter dependence of the ordinary bare loop expansion 
by the same resummation in HTL propagators and vertices, when calculating damping rates and plasma 
frequencies, represents an important check on its physical consistency [9, 14, 17-21]. 

The technical evaluation of the two-loop diagrams in the imaginary time formalism [3, 5] is facilitated 
considerably by using Gaudin's method [22, 23], which we review in Appendix A. This very efficient 
method of performing the frequency sums makes calculations of diagrams involving many finite temperature 
propagators feasible. The Gaudin method also makes the hard-soft factorization pattern at two-loops more 
evident, and makes feasible the analysis of dressed and/or even higher loop diagrams, at least for the 
purpose of extracting their leading order behavior in the high temperature limit. 

The outline of the paper is as follows. In Section II we review the one- loop finite temperature electron 
self-energy, including both the LO HTL and NLO terms, and the latter's dependence upon the gauge fixing 
parameter £. In Section III we present the basic calculation of the three diagrams of two different topologies 
in the two-loop electron self-energy by the Gaudin method in the Feynman gauge, introducing the photon 
mass regulator for the most infrared sensitive bubble diagram. In Section IV we present the detailed 
analysis of the two-loop integrals obtained in Section III, for a fermion at rest q = 0. We show that the 
leading contributions in the high temperature expansion of the two-loop diagrams of both topologies come 
from one loop hard, the other soft, thus demonstrating explicitly that there are no hard-hard contributions, 
and the soft-soft contributions are subleading at this order. In Section V we compare our results with the 
insertions of effective HTL self-energies and vertices in the one- loop S, explicitly verifying the equivalence 
of these insertions for both the real and imaginary parts to order e 2 T, which is necessary for the HTL 
resummation program to work. In Section VI we study the gauge parameter dependence of the two- 
loopself-energy, and show that it does not cancel even on shell, in neither its real or imaginary parts. 
Section VII contains a summary of the main results for the two-loop self-energy to order e 2 T, concluding 
with a discussion of the lessons for a systematic reorganization of the loop expansion in high temperature 
gauge theories to higher orders. 
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FIG. 1: One- loop contribution to the fermion self-energy. 

II. THE ONE-LOOP ELECTRON SELF-ENERGY 

We begin by setting our conventions and reviewing the existing leading order (LO) [3, 24] and NLO 
[10, 11] calculations of the one-loop QED electron finite temperature self-energy, 



Si(Q) = -ie 2 [ llx G{Q - K) lv D^{K) (2.1) 



represented in Fig. 1. For the purposes of this paper, we use the Feynman rule conventions of [25] wherein 
the tree-level electron-photon vertex is —iej^. In (2.1) the tree-level fermion propagator is given by 



G(K) = -1 = fiD(K) , D(K) = i/(K 2 + ie) , (2.2) 

with f[ = j^Ku = &o7o — k • 7. The tree-level photon propagator in the one-parameter family of covariant 
^-gauges is 



(2.3) 



We employ the somewhat redundant notation of D(K) and D(K) for the same function of K 2 in order 
to distinguish propagators with fermionic Matsubara frequencies from those with bosonic Matsubara fre- 
quencies in the imaginary time formalism of finite temperature field theory [3, 5]. Thus, in imaginary time 
where the generic four-momentum becomes = (iu} n ,\c), 

— ^ 

D(iu n , k) = -iA(iu n , k) = — — , uj n = 2irnT ; (2.4a) 

\lU)n ) T K 

D(iu n , k) = -i£(w w , k) = ~* w n = (2n + 1)ttT , (2.4b) 

The indices of the four-momentum are raised and lowered with the metric tensor g^ v = diag(l, —1, —1, — 1). 
The Dirac 7-matrices satisfy the anticommutation relation {7^, 7^} = 2g^ u , so that (7 )' = 7 is Hermitian. 
For evaluation of loop integrals such as that in (2.1) we use the convention 
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Moreover, our conventions are such that the electron and photon proper self-energies, £ and 11^, contribute 
to the full inverse propagators in the form 

ig- 1 (Q) = iG-\Q)-Z(Q), iV-^Q) = iD^(Q)-U^(Q), (2.6) 

for the electron and photon respectively. The sign convention for £ here is the opposite of that of [11]. 
It will also be convenient to work with the spectral representation of the propagators, which is given by 

D^(iu n ,k) = -i d ^Pph3. = _,AH-n,k), (2.7a) 
J-oo 2vr fc - iu n 

G(iu n ,k) = -i r d ^PJ^l = -iA F (iu n ,V, (2.7b) 
J-oo 2vr k - iuj n 

where 

p^(k ,k) = - (V + (1 - OK^-^j p(k ,k), (2.8a) 

PF (k ,k) = fip(k ,k). (2.8b) 

Here p(ko, k) is the spectral function of a free massless boson, 

p{K) = p(k , k) = 2ire(k )5(k 2 - k 2 ) , (2.9) 

and e(ko) = 9(ko) — 9(—ko) is the sign function. Accordingly, the spectral representation of A(io; n ,k) is 

A(lU n ,k) = - — 2=/ — : , 2.10 

K 2 J-oo 2vr k - iu} n 

and similarly for A(iu n , k), with the only difference that in this latter case oj n is fermionic. 

A. Feynman Gauge Evaluation 

We consider first the Feynman gauge condition (£ = 1). In this case, using the spectral representation 
of the propagators, the self-energy (2.1) for a fermion with external Matsubara frequency v q becomes 



ExK, q) = -2e 2 SJl'tjl T^P(k , k)p F ( Po , q - k) T £ 



1 



1 



_P0 ~i(Vq- "r. 



(2.11) 



The frequency sum can be performed using Gaudin's method [22], reviewed recently in [23] and in Ap- 
pendix A. 

Gaudin's method relies first on the spectral representation of finite temperature propagators in the 
imaginary time formalism (2.10), and second, on a graphical representation of a particular partial fraction- 
ing decomposition of the products of denominators that arise in finite temperature perturbation theory into 
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a sum of products, in each of which the only Matsubara frequencies which appear are the independent ones 
to be summed. Since an £-\oop finite temperature diagram contains £ independent Matsubara sums, this 
means that it can be written as a sum of terms each one of which contains exactly t independent frequency 
sums, which can be performed simply, no matter how many propagators or vertices the original Feynman 
diagram contains. This disentangling of the independent and dependent Matsubara frequencies greatly 
facilitates the computation of multi-loop finite temperature diagrams, or even one-loop finite temperature 
diagrams with multiple vertex insertions. 

For the one-loop diagram of Fig. 1 the Gaudin decomposition into terms with the one independent 
Matsubara frequency to be summed amounts to the simple partial fraction identity, 



k - iuj n 

e iVqT2 

ko+PO- Wq 



e i(u q -Ul n )T 2 

E jw n (ri-r 2 ) 
: 
m k - tu n 



+ 



ko+Po- iv q 



e «(^-w n )(T 2 -ri) 
PO ~ i{Vq ~ W n ) 



(2.12) 



so that only one denominator in each term now contains the Matsubara frequency to be summed. Each 
of the Tk are infinitesimals, assumed positive, assigned to each line; for example Tk = ke with e — > + , so 
that Ti — T2 has a definite sign. Shifting then the summation over the bosonic frequency ui n = 27rnT in the 
second term of (2.12) to the fermionic frequency Cj n = v q — uj n = (2n+ 1)ttT, and applying the fundamental 
summation formulas (A3) and (A4) for bosonic and fermionic frequencies respectively, we obtain for the 
sum in (2.11) or (2.12) the result 

1 n F (p ) + ns(-fco) 



(2.13) 



(fc - iUn) [PO -i{ v q- Un)] k + p ~ W q 

in the limit — > + with n — T2 — > 0~. Here nB(E),np(E) are the Bose-Einstein and Fermi-Dirac 
distribution functions given by 

1 1 



n B (E) = 



n F (E) = 



(2.14) 



exp(E/T) - 1 ' v ' exp(£/T) + 1 

The power of the Gaudin method lies in the fact that once the propagator denominators with independent 
Matsubara frequencies are isolated, as in (2.12), one never needs any other Matsubara frequency sums 
other than the fundamental ones (A3) and (A4), no matter how many propagators or vertices the original 
Feynman diagram contains. This advantage becomes more and more evident at higher loop orders, or in 
diagrams with larger numbers of finite temperature propagators. 

We remark also that had we assigned the r» differently, so that, e.g. t\ — ti — >■ + , we would have 
obtained +[np(— pq) + ne(/co)] instead of — [wf(po) + n B(— ko)} in the numerator of (2.13). However these 
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two expressions are equal, by use of the identities 

n B (-E) = -l-n B (E), n F {-E) = 1 - n F (E), (2.15) 

for the Bose-Einstein and Fermi-Dirac statistical distributions. Hence although one needs some regulator 
to define the frequency sums (A3) and (A4), appearing in intermediate steps in the Gaudin method, one 
can check that the results for any diagram containing at least two Matsubara propagators (whose frequency 
sums converge absolutely) are in fact independent of the regulator used. 

Performing next the frequency integrals over ko and po in (2.11), by making use of the 5-functions in 
the spectral representations (2.8) and (2.9), and using the first identity in (2.15), we obtain 



EiK,q)=2e 2 /t^T^ £ [(t°^i - (q - k) • 7 Sl s 2 ) 

s 1 ,s 2 =±l 



1 + n B (siEi) - n F (s 2 E 2 ) 

lUn — S\E\ — S2E2 



, (2.16) 



where E\ = k, E 2 = |q — k|. For each particular value of s± and s 2 the terms in the sum can be associated 
with physical processes like pair creation and annihilation or scattering off of particles in the heat bath 
[26]. 

Examining (2.16) one sees that for soft external four-momentum there are two types of terms: one with 
s 2 = —si and the other with s 2 = s\. In the first case the energy denominator is small for |q| <C |k|, 
since E\ — E 2 c± k ■ q/fc, and the k integration is dominated by hard k = |k| ~ T. These hard momentum 
integrals are proportional to T 2 and give the LO result in the HTL approximation. Physically these 
terms correspond to Landau damping in the plasma, which includes processes possible only at nonzero 
temperature: particle absorption from the heat bath and emission into the heat bath. 

In the second case, that is with s 2 = s\ = ±1, the two energies E\ and E 2 add rather than cancel 
for large k in the denominator of (2.16), and hence there are additional powers of k in the denominator 
which makes the integral over k more sensitive to its small k range. Since 1 + nB(s\E\) — n F (s\E 2 ) = 
e{si)(l + nB{E\) — n F (E 2 )), the appearance of the "1" shows that here one encounters processes which are 
possible also at zero temperature, e.g. pair creation. The NLO terms in the high temperature expansion 
of the one-loop self-energy come from these types of processes. 

Focusing on a fermion with vanishing spatial three-momentum (q = 0) one has E\ = E 2 = k and the 
above two possibilities give 

LO HTL: from V —>-—\n B (k)+n F (k)]—, (2.17) 



2k 1 Hv q 

S2 = — S\=±l H 



NLO HTL: from £ ->-*. [l + „„(*) - »,(*)] (— + — L^) . 

S2=Sl=±l 



(2.18) 
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Using (2.17) and (2.18) one obtains then 

e 2 f°° e 2 iu f°° k 

Si(^ g ,q = 0) = ^tPw 70 J dkk [ n B( k ) + n F(k)} - -^-Jo J dk ^ 2 - — [ l + n B {k) - n F (k)] . 

(2-19) 

In the last expression the "1" in square brackets is the logarithmically UV divergent vacuum contribution, 
which must be renormalized in the standard manner. Since we are interested only in the high T dependence 
of the self-energy in the hot plasma, we may discard such vacuum contributions. It is important to point 
out that because of the identities (2.15), identifying the vacuum contribution unambiguously requires first 
expressing all statistical distributions as functions of energy arguments which are positive. This can be 
done most conveniently before doing the frequency integrals by using the following expressions: 

n B (E) = -6(-E) + e(E)n B (\E\), n B {-E) = -6(E) - e(E)n B (\E\), 

(2.20) 

n F (E) = 6(—E) + e(E)n F (\E\), n F (-E) = 9(E) - e(E)n F (\E\). 

Clearly the k integral in the first term of (2.19) is hard in the sense that a positive power of k in the 
integrand is cut off only by a statistical distribution function n B or n F at k ~ T, and we obtain 

2 T 2 2 T 2 2 T 2 

dkk[n B (k)+n F (k)] =— + — = —. (2.21) 

In contrast, in the second integral of (2.19) the integrand is a decreasing function of k for large k, even 
without the statistical distribution functions. Integrals of this kind are soft, in the sense that the small 
k range is important, and the n B or n F thermal functions serve only to cut off at most an otherwise 
logarithmically diverging integral at large k, in which the argument of the logarithm is sensitive to the soft 
momentum region. The high temperature limit of the second soft term in (2.19) is evaluated in Appendix 
B, where it is shown that 



f°° k 

dk-^ — - — ry— n F (k) 



l lu( T ) • .... (2.22a) 



J 

JO 



o p - {iv q y/A 



+ -)+... , (2.22b) 



= W +iO+ 2 \ u 



where only the high temperature asymptotic values of the integrals are displayed, neglecting terms that 
do not grow with T. Thus, in Feynman gauge, 

e 2 T 2 e 2 ( T \ 
£i(w,0)~— 10 + — 10 ( u ] n --i ir T) , (£ = 1). (2.23) 



The first term is the leading order (LO) HTL result, which receives contributions from both the electron and 
photon thermal distributions in (2.21). It is purely real, and in fact, independent of the gauge parameter 
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£. The second two terms in (2.23) are the next-to-leading order (NLO) or subleading terms in the high 
temperature expansion, which come from the soft momentum region in the integrals (2.22), and yield both 
a real and imaginary part. The contribution of the Fermi-Dirac statistical factor to the imaginary part and 
of the zero temperature vacuum part are subleading to this (NNLO) in the high temperature expansion. 

The real part of the off-shell one-loop self-energy Si in (2.23) illustrates the infrared divergence at 
small external frequencies co — > 0. The on-shell condition at LO is obtained by finding the zeroes of the 
one-loop corrected inverse propagator iG~ l (Q) — in (2.6). For vanishing spatial momentum q = 0, 

this on-shell condition is 

w 70 = S(w,q = 0) (2.24) 

or to leading order, 



e 



T 2 



"--5T + - < 2 ' 25 > 

where the ellipsis denotes the NLO terms. Thus on the fermion mass shell the tree level inverse Dirac 
propagator is set equal to the one-loop result and uj ~ eT/y/8 [3, 27, 28]. It is the fact that oj ~ eT on-shell 
owing to the infrared enhancement of the one-loop self-energy at small oj that leads to inverse powers of 
the coupling at higher loop orders, re-ordering of the perturbative expansion, and eventually the need for 
further resummation. 



B. Gauge Parameter Dependence of Si 

The general expression for the one-loop self-energy at NLO in the HTL approximation and for q 7^ 0, in 
a general £ 7^ 1 gauge can be found in [10] (real part) and [11] (both real and imaginary parts). Explicitly, 
the £ dependent term is 

f 1 K fM K u 

^i(Q) = (l-flie^—^ — (2.26) 



Making the substitution = — fi) + $ we can rewrite this as 

" . 1 



^ l( g) = (i-o^/^ 



<&-{<& + *) 



(2.27) 



{K 2 ) 2 

Since the very last term here involving only ]/C vanishes by symmetry, in this form one can see that if 
the lowest order tree level on-shell condition for a massless fermion $ = is used, the gauge parameter 
dependence of the one-loop self-energy vanishes, provided that there are no infrared divergences in this 
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limit. This is indeed the result that one obtains in the well-known QED perturbation theory at zero 
temperature with a finite electron mass. 

At high temperatures, the existence of the w _1 infrared behavior of the LO HTL self-energy in (2.23), 
and imposition of the modified HTL on-shell condition (2.25) implies that ($ = lu^/o / for a fermion at 
rest in the plasma. Hence (2.27) no longer vanishes for HTL dressed fermions and the gauge parameter 
independence of the modified fermion dispersion relation at any given loop order is no longer guaranteed, 
even on-shell. This is clearly the result of using an on-shell condition such as (2.25) which mixes loop 
orders. 

To evaluate ^Si(Q) explicitly in the high temperature limit, we group the terms in a somewhat different 
way. Using the fact that fC($- fC)fi = Q 2 fi- K 2 $-(Q- K) 2 fi, and that ^ J K fiD M 2 (K) vanishes, 
we may also write ^Si(Q) in the form 

^ 1 (Q) = (l-Oie 2 {-iQ 2 ^ J K fiD M2 (K)D(Q-K) ^ + i® J^D(K)D(Q - #)} , (2.28) 

where 

D m*W = k2=m2 ( 2 - 29 ) 
is the bosonic propagator with finite M 2 . The corresponding spectral representation is 

n (K\ „■ f°° dk p M 2(k ,k) 

D M 2{K) = -i : , (2.30) 

J-oo 27r h - lUJ n 

with 

p M2 (fc ,k) = 2ire(k )6(k% - k 2 - M 2 ) . (2.31) 

The technique of introducing a photon mass M, differentiating with respect to it, and taking the M — > 
limit at the end of the calculation is a useful method for defining the double pole term, (K 2 )~ 2 appearing 
in both ^Ei and the two- loop self-energy bubble diagram of the next section. At this point the photon 
mass M is introduced as a mathematical convenience and infrared regulator only, and the limit M — > 
should be taken at the end of the calculation. A finite photon self-energy of order M ~ eT is obtained 
only after HTL resummation. 

The reason for the rearrangement of ^Si into the form (2.28) is that each of the terms in (2.28) are 
finite in the limit M — > 0, whereas as we shall show in Sec. VI, each of the terms in (2.27) is M dependent 
and divergent in this limit, with the sum of terms in ^Si finite. In imaginary time, using the spectral 
representation of the propagators (2.10) and (2.30) we have 



*£i(u/g,q) = (l-£)e 2 7o x 
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M 2 =0 



\ - Q W LL ho PM2 ( "°' k) L to p ^ q - k)T e (A.-^bo-^,-^)] 

+0 / r ? r ^o,k) P ( Po ,q-k)T^ , 1 J. (2.32) 

Since ui n is bosonic and v q —oj n is fermionic, the Matsubara sum here is identical to that given in (2.13). One 
next performs the frequency integrals over ko and po in (2.32). For the integrals over the frequencies which 
are not arguments of the statistical factors, the spectral decomposition (2.10) may be undone, re-obtaining 
Matsubara propagators. The remaining frequency integrals are performed using the Dirac 5 functions 
in the spectral functions after using the relations (2.20), which allows the separation of the vacuum and 
thermal pieces by making the arguments of the statistical distributions explicitly positive. Keeping only 
the thermal part and taking q = 0, in this way one obtains 

(i-Oe 2 ., 2 d [2 f°° dk kn F (k) 2 r°° d] k 2 E k n B (E k ) 



^K,q = 0) = ^2—70 K) ^ J o 




^ 2 -/i 2 iv q Jo ' E%-»% j M2=Q 
^■(■^ M^-Mfc)]}, (2-33) 

where E\ = k 2 + M 2 and \i± = {{iv q ) 2 ± M 2 )/(2iu q ). We have omitted from the square bracket terms 
independent of M 2 which give a vanishing contribution upon differentiation. 

The last integral in (2.33) is the same one obtained in the Feynman gauge calculation, c.f. (2.19). After 
analytical continuation (iv q — > oo + i0 + ), the ^-dependent real part of the self-energy comes only from this 
term. From (2.22) it also contributes to the imaginary part. The first integral of (2.33) is of the form given 
in (Bl) with the replacement ia — > H-/T. Using (BlOa) one can readily check that the term containing this 
integral will give a contribution which is subleading to (2.34). The second integral in (2.33) contributes to 
the imaginary part of ^Si as shown by (B24). All these £ dependent terms arise from the soft momentum 
region. The final result for the high temperature expansion of the £ dependent part of the one-loop fermion 
self-energy is 

€£1(07,0) = (1 -0^70 (""In - + —) ■ (2.34) 

Unlike the LO HTL result, the NLO terms and therefore the imaginary part of the one-loop self-energy for 
a fermion at rest in the plasma is gauge-parameter dependent, even after the LO HTL on-shell condition 
(2.25) has been utilized. This clearly indicates that the imaginary part obtained from the NLO one-loop 
self-energy calculation alone cannot be the full answer for the dressed quasi-particle damping rate, and 
both the real and imaginary parts of S should receive additional contributions of the same magnitude in 
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Q-P,5 



P,5 



v 





Q-K,2 



Q-PA 



Q-P,3 



P + K,3 



Q-P-K,3 



FIG. 2: Two-loop contributions to the fermion self-energy: from left to right, the bubble (B) diagram, rainbow (R) 
diagram, and the crossed photon (C) diagram. 

the weak coupling e at higher loops, in order to obtain a physical result independent of £. We show next 
that terms of the same or larger magnitude to (2.34) are indeed generated at two-loop order, for to ~ eT 
due to increased sensitivity to the infrared of the two-loop self-energy. 



In this section we compute the two-loop electron self-energy in the Feynman gauge £ = 1, perform 
the Matsubara sums in the imaginary time formalism by Gaudin's method, and reduce the expressions 
to integrals over two spatial loop momenta k and p, for general external electron energy and momentum 
(iv q ,q). The leading terms in the high temperature expansion arise from only those terms with two 
statistical factors, ns or np. The analysis of the integrals, the explicit extraction of their leading high 
temperature, most infrared sensitive terms, and the analytic continuation to real electron energies iv q — > 
to + i0 + for q = electrons at rest in the plasma are performed in the next section. 

At two-loops the electron self-energy is composed of the three diagrams pictured in Fig. 2, which we call 
the bubble (B) diagram, the rainbow (R) diagram, and the crossed photon (C) diagram. The corresponding 
contributions to £ will be denoted by S B , S R and S c respectively. 

An earlier partial calculation of the two-loop electron self-energy at finite temperatures has been given 
the literature in Ref. [12]. The authors of this work were interested in the ultraviolet renormalizability 
of QED and electron mass renormalization at finite temperatures. For that reason they used real time 
Feynman propagators (rather than the appropriate causal real time Schwinger-Keldysh formalism), and 
did not calculate the imaginary parts of any of the two- loop diagrams. They also did not calculate the 
contribution of the bubble diagram (B), which turns out to give the dominant contribution to the real part 
at leading order in the high temperature expansion, in which we are interested. 

The calculation of such multi-loop diagrams in the imaginary time formalism is greatly facilitated 
by Gaudin's method for performing the Matsubara sums [22]. Since each of the two-loop self-energy 



III. THE TWO-LOOP ELECTRON SELF-ENERGY 
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diagrams contains five finite temperature propagators, finding the appropriate partial fraction identity 
analogous to (2.12) for disentangling the two independent Matsubara frequencies is already quite non- 
trivial. Gaudin's method for doing this involves a graphical decomposition of each original two-loop 
Feynman diagram in E into a sum of tree diagrams, each containing three propagators, and the complement 
of the tree containing the other two propagators which carry the two independent Matsubara frequencies. 
This graphical decomposition is illustrated in Fig. 3, for the case of the bubble (B) self-energy diagram. 
Calculations in scalar field theory using Gaudin's method were presented in [23] on a two-loop diagram 
having the topology of diagram (C). Diagrams of type (R) or (B) were outside the interest of the authors 
of [23] because, having a self-energy insertion on an internal line, they are not in the class of two-particle 
irreducible (2PI) diagrams. The method works equally well on diagrams of all topologies. However, 
since two fermion propagators (diagram R) or two photon propagators (diagram B) have the same four- 
momentum, the two diagrams suffer from a double pole problem (in any method). In Gaudin's method 
this problem manifest itself in two different ways, as one sees from the Gaudin tree graphs of Fig. 3. The 
first way is illustrated by the first four graphs, where one encounters the product of a spectral function 
corresponding to a line which acquires a statistical factor and a Matsubara propagator with the same four- 
momentum, that is 5(pq — p 2 ) / (p 2 — p 2 ,). The second way in which the double pole problem arises is found 
in the last three Gaudin tree graphs, where one encounters the product of two Matsubara propagators 
sharing the same 4-momentum. In both cases the expressions can be regularized by a massive regulator, 
which effectively defines ill-defined double pole quantities by the replacements 

p(p°,p)A(p»,p)^ ^PmV,p)= r 2 ^{^(P°)s((P y 2 - P 2 - M*)}, (3.1a) 

The M 2 ->• limit can be examined then at the end of the calculation. As a technical remark, M 
may be considered to have an infinitesimal negative imaginary part M 2 — > M 2 — i0 + , if necessary in any 
intermediate steps of the calculation. 

A. The Bubble Diagram 

The contribution of the bubble graph of Fig. 2 reads 

S B (Q) = ~ie A [ [ 7 aG(Q - P) lp D^{P)ti [ 7f4 G(P + K) lv G{K)\ D^(P). (3.2) 

JK JP 
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FIG. 3: The Gaudin tree graphs corresponding to the bubble diagram of Fig. 2. The thin lines carry the independent 
Matsubara frequencies while the thick lines belong to the tree (see Appendix A) . The assignment of the momenta 
of the lines are the same as in the corresponding bubble Feynman diagram of Fig. 2. 

The double pole of this diagram leads both to possible infrared and collinear divergences (see e.g. [29] 
and references therein), which requires regularization. For this reason we introduce a mass M for the 
photon which will be kept to the end of the calculation, in order both to regulate the divergences of the 
perturbative two-loopself-energy and also to estimate the coupling and temperature dependence of the 
resummed electron self-energy when M is replaced by the Debye screening mass ~ eT. 
In imaginary time, using the spectral representation of the propagators, we have 

5 - '„o\ 



£ B ("Vq) 




J a pF(P 5 h(3P°""(Pi)p^(P4)tT{j ll p F (Psh u p F (P2)} 



Z^ll o 

n,m i=l * 1 



lUJi 



(3-3) 



where Pi = P 4 = P, P 2 = K, P 3 = P + K, and P 5 = Q - P with Q = (iv q ,q), P = (ioJ n ,p), and 
K = (iu; m ,k). The Matsubara sums performed with Gaudin's method gives for the last factor of (3.3), 



t 2 VTT 



i 



n,m t=l 1 



ILUi 



2 E 

n,m 



1 



pi - iu. 



1 



P2 - iuJ m 

nF{vV)nB{-vlV 

P2+P4- P3 



1 



P3 



lOJr, 



(Pi -PD(P$ +Pl - iv q) 

1 

(Pi -P°)(P° +Pl - iy q) 

n F {pl)n F {pl) 



+ 



lUJr, 



1 



1 



1 



PI - iuj n J \pl - iv q + iuj n 
nB{-Pi)n F (pl)^ 



n F (p2)n B (pT) 
P3-P2- Pa 



{pI-Pi){pI+Pi ~ iv q ) 
1 

(P° -Pl)($>l +Pl ~ iv q) 



P2+P1- P3 
nB{Pi)n F {p Q 2 ) 

pI~p\-pI J 



+ 



+ 



(P? + P°2 ~ P3) (P° + P2 ~ P3)(P°5 +P3-P2- ™<1 
1 



(P? + P5 - iv q)iPA +P5- iu q)(P2 -P3-P5 + 



n F {p^)n F (-p%) + n F (pl)n F {pl) 



q> 



(3.4) 



where each of the seven terms corresponds respectively to each of the seven Gaudin tree graphs shown in 
Fig. 3. 
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For each of these terms two of the lines (the thin lines in Fig. 3) have associated statistical factors n B 
or n F , while the remaining three lines (the thick lines in Fig. 3) have no statistical factors associated with 
them. The integrations over the p® for these latter three propagators can be performed simply by undoing 
the substitution of the spectral representation (2.7a) and (2.7b) in favor of the original imaginary time 
propagators, A^ v or A F . Then the first and second of the seven terms give identical contributions with p\ 
replaced by p®, as do the third and fourth terms with the same replacement. Dropping the superscript 
on the remaining integration variables p\ to simplify the notation, we obtain then from (3.3) and (3.4) 

S B (^g,q) = e 4 ^y ^2 J ^ J d ^n F (p z )n B {-p A )x 

-y a A F (iv q - p 4 , q - p)7 /3 A^(p 4 , p)p^(p4, p)tr {i^Pf(P3, P + k) 7l ,A F (p 3 - p 4 , k)} 

~j a A F {iv q - p 4 , q - p)7 /3 A^(p 4 , p)^(p 4 , p)tr {j^A F (p 2 + p 4 , P + k)j u p F (p 2 , k)} 

j a A F (iis q +p 2 -J?3,q- p)7,3A aM (p3 - p 2 ,p)A" /3 (p 3 -p 2 ,p)tr {7mPf(P3, P + k)7!,pF(P2,k)} 

dp 2 /" dp 5 



7aPF(P5,q- p)7/?A a ^(ii/g -p 5 ,p)A^(ii/ g -p 5 ,p)tr{7 (U A jF (ii/ 9 +p 2 - f>5, P + ^)1vPf{P2, k)} 

lapF{Pb, q - p)7/3A^(ii/g - p 5 , p)A^(wg - p 5 , p)tr {7 m pf(P3, P + k)7„A F (p 3 + p 5 - iv q , k)}} .(3.5) 

Expressing the fermionic and gauge Euclidean propagators A^ and A Miy , and the spectral densities pp and 
p fiV in terms of the bosonic forms A and p defined in (2.10) and (2.9), and evaluating the expression (3.5) 
in the £ = 1 gauge, we obtain the covariant form for the bubble diagram, 

E.K, q)U = ^/ t / p £f£f^WP)x 

{2 [n F (-k°)n B (-p°) + n F (k )n B (p°j\ B(K; K + P-Q - P)A(P)A(K + P)A(Q - P) 

+ n F (k )n F (p°)B(K; P;K + Q- P)A 2 (K - P)A{K + Q - P) 
- n F (k°)np(-p°) + n F (-k°)n F (p°)j B(K; K + Q - P; P)A 2 {Q - P)A(K + Q- P)J , (3.6) 

where we have defined the factor arising from the Dirac algebra, 

B(P 2 ;P 3 ;P 5 ) ee B(p° 2 ,p 2 ;plp 3 ;plp 5 ) ee 57 "tr { 7M f 3 7,f 2 } = 8[f 2 (P 3 ■ ft,) + f 3 (P 2 ■ ft)] , (3.7) 

and relabeled in the first term of (3.5) p 3 — > —k°,p4 — > p°,~k — > — k — p ; in the second term of (3.5) 
P2 — > k°,p4: — > p°; in the third term of (3.5) p 2 — > k°,p 3 — > p°,p — > p — k; in the fourth term of (3.5) 
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P2 — > A: ,p5 — > p°, p — > — p + q; and finally in the fifth term of (3.5), ps — > — k ,ps — > p°, p — > — p + q, k — > 
— k + p — q, in order to arrive at (3.6). 

From the form of (3.6) we observe the double pole structure of the bubble diagram as discussed at 
the beginning of this section. The first term of (3.6) within curly brackets contains the singular product 
p(P)A(P), which is regulated by the replacement (3.1a). This term will turn out to give the most singular 
contribution to the twodoop S in this limit. The second and third terms of (3.6) within the curly brackets 
involve the square of a zero temperature bosonic propagator A 2 (K — P) or A 2 (Q — P), both of which are 
regulated by the replacement (3.1b). 

In order to perform the detailed analysis of the integrals appearing (3.6) and exhibit explicitly their 
hard-soft pattern, it turns out to be more convenient to express (3.6) in a slightly different form. Leaving 
the first term as it is, in the second term of (3.6) we shift P — > P+K (i.e. p° — > p°+k° and p — > p+k), while 
in the third term of (3.6) we undo the shifts defined in the passage from (3.5) to (3.6) which were necessary 
in order to express (3.5) in the compact covariant form (3.6), that is, simply relabel p2 = k ,p$ = p° and 
P3 = /c°,p5 = p° in the last two terms of (3.5), without changing the integration variables p and k. After 
introducing the M 2 regulator (3.1a) of the double pole terms, and evaluating the k° integrals, one obtains 
from either (3.5) or (3.6) the three terms 



E B (iu q ,q)\ 



T, B1 (iu q ,q) + £ B2 (iz^,q; 



.) + Y, B3 (iv q ,q), 



(3.8) 



with 




oo 



dp 1 
~2ir2k 



£ B i(^ ? ,q) 



4 v u 
e hm — — j 

m 2 ^o dM 2 



k J P 



A(iu q -p ,q-p)p M 2(p°,p) r< A(r/c+p°,k + p) 



oo 



r=±l 



T, B2 (iv q ,q) 




n F (rk)n F (p° + rk)B(rk, k; p° + rk, p + k; iu q - p°, q - p) \ , (3.9b) 



£ B3 (iz/ ? ,q) 



Q r r r-CG ^ ( 

e4 J^o am L I L ^ p(p0 ' q " p) Am2 {Wq - p) 5/1 




nF(-p°)B(rk, k;iv q + rk - p°, p + k;p°, q - p)A(iv q + rk —p°,p + k) 
P ,t^ n F {p Q )B{rk +p°- iu q , k; rk, p + k; p°, q - p)A{rk + p° - iv q , k) \ , (3.9c) 
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where k = \k\,p = |p|. We will evaluate the leading order terms in the high temperature expansion of (3.9) 
for q = in the next section. 



B. The Rainbow Diagram 

The contribution of the rainbow (R) graph of Fig. 2 is given by 

£ R (Q)=W [ 7 a G(P 1 )^G(P 3 ) 7u D^(P 2 )G(P^D^(P 5 ), (3.10) 

JK JP 

where we have introduced the labels, Pi for the five different four-momenta, P\ = P4 = Q — P, P 2 = K, 
P3 = Q — P — K, and P5 = P, with Q = (iu q ,q), P = (iu) n ,p), and K = (iw m ,k), as in the previous B 
diagram. 

It is possible to evaluate this rainbow diagram in a completely analogous manner as the bubble diagram. 
Since the topology of the two graphs is the same, one will obtain seven tree diagrams by the Gaudin method 
for the rainbow self-energy, analogous to the seven tree diagrams of Fig. 3. In the case of the rainbow 
diagram the differing Dirac structure allows for some simplification. In the Feynman gauge £ = 1, we have 

7af i7„f 3 7 M f 4 7 a = 4f 4 f 3 f l = 4 { [Pi + Pi - (P 1 - P 3 ) 2 ] f 4 - Plf 3} - (3-H) 



since P\ = P4. Using also that Pi — P3 = P2 we have 
7af i7^ 3 7 M f 4 7 a _ 



p2 p2 p2 p2 p2 p2 p2 p2 
L^2 ^3 ^4 ^5 ^2 ^3 ^4 ^5 J 



+4 (il4G?-i>' (3 - i2) 



As in the bubble diagram there is a potential double pole problem which shows up only in the second set 
of terms in (3.12), and can most conveniently handled by introducing a fermion mass m, and making the 
replacement 

»f 4 9 ( if A \_ dG m *(P 4 ) 

(P|)2 ^ dm 2 VP 4 2 -m 2 J~ dm 2 ' 1 J 

and taking the limit m 2 — > at the end of the calculation. Because this is a fermionic double pole rather 
than the bosonic one in the bubble diagram, the limit will be less singular, and in fact not give rise to 
divergences in the leading order of the high temperature expansion. 

The two sets of terms in (3.12), denoted by Rl and R2 give for the full rainbow two- loop self-energy 
the expression 

S R (Q)| ?=1 = S R1 (Q) + S R2 (Q) , (3.14) 
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Q-P,l Q-P-K,3 

FIG. 4: The diagrammatic representation of the integrals in (3.15a) obtained by reducing the original rainbow diagram 
and its decomposition in Gaudin tree graphs. The dashed (solid) line means bosonic (fermionic) propagator. As 
in Fig. 3 the two thin lines carry the two independent Matsubara frequencies to be summed, while the thick lines 
belong to the Gaudin tree and represent propagators that contain the time component variable of the appropriate 
momenta p° in place of the Matsubara frequencies. 



with 



Ski(Q) 



-4e q 



K JP 



£ R2 (Q) = 4ie 4 ( I 



D(P 2 )D(P 3 )G(P 4 )D(P 5 ) - D(P 2 )G(P 3 )D(P 4 )D(P 5 ) 
d 



K D(K) ) a^j P G ^ p)D ^- p) 



m 2 =0 



(3.15a) 
(3.15b) 



To obtain the second R2 term we have used for the last term inside the parentheses of (3.12) the shift 
K —7- Q — P — K and then for both terms within the parentheses the shift P — >• Q — P. The advantage 
of the form (3.15) is that whereas the original expression (3.10) has five propagators, the Rl terms have 
only four propagators, and the R2 terms factorize into the product of a difference of two tadpoles and a 
derivative of a one-loop self-energy with only two propagators. In these expressions it is important to keep 
track of which momenta are bosonic and which are fermionic. The latter propagators are denoted with a 
tilde. 

In imaginary time, using the spectral representation of the propagators, we have for the two contribu- 
tions 



£ R1 (^,q)=4e 4 / / J] 
JkJp \ i=2 



2vr 



[p(P3, P3)pf(pI, P4) - Pf(pI, P3)p(pI, P4) 



p(pg,P2Mpg,P5) r 2 En 



— lUJj 
n,mi=2^ , 



(3.16a) 



£ R2 (^ g ,q 




4e 4 T 



dm 2 



^ k \ m m ) 

f f°° dp f°° dk° , n r, w,n x^v- 1 



.(3.16b) 



m 2 =0 



The Matsubara frequency sums in the Rl terms may now be performed by Gaudin's method. The Gaudin 
relevant tree diagrams represented diagrammatically in Fig. 4. Having one propagator less compared to 
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the original graph, it is somewhat easier to compute the sums, and we obtain 



n B {-p° 2 )n F (-p° 3 ) n F {-pl)n B {-p° 2 ) 



, , ,„ , =2 P°i ~ iu>i (p° 4 -P° 2 - P° 3 )(p° 5 +P° 2 +P° 3 - iv q ) (p§ - p\ + pg) (pg + p2 " w?) 

n F (-^)n F (P3) "-B(-P2) n B(P5) 



(P°2 -P°A+ Ps) (P°5 +P4- Wq) (P° + ?5 " (P° +^2+^5" 



(3.17) 



(P° +P°5~ iv l)(P2 +P3+P5- iv l) ' 
with each term corresponding to a Gaudin tree graph of Fig. 4. In (3.16b), u n is fermionic and v q — co n is 

bosonic, and we obtain 

T V- 1 = n F (-p°) + n B (r°) 

^ (p° - iuo n )(r - i(u q - u n )) p° + r°-iu q 

for the Matsubara sum. 

Substituting (3.17) into (3.16a), the integrations over the two p\ which do not appear in the statistical 
distributions can be performed simply by using the definition (2.10) of the spectral representation of the 
Euclidean propagator function A. By relabeling p\ — > k°,p® — > p°, p — > — p + q — k in the first of the five 
terms arising from (3.17), p\ — > k°,p\ — > p°, p — > — p + q in the second term, p® — > k Q ,p\ — > p°, k — > — k + 
p,p —7- — p + qin the third term, p 2 — > k°,p® — > p in the fourth term, and p® — > k°,p® — > p°, k — > — k+q— p 
in the fifth term arising from (3.17), we arrive at the covariant form, 

S R1 (^,q) = -4e 4 ^| J°° ^ 1°° ^^n B (-k )n F (-p°)[A(K + P)A(Q - K - P) 

+ A(P - K)A(Q - P)] + (f - fi)n F {k°)n F (-p°)A{P - K)A{Q - P) 

+ n B (p°)A(Q - P)A(Q -P-K) [f n B (-k°) + (# - f - fi)n F (k )] |p(fc°, k)p(p°, p) , (3.19) 

for (3.16a). 

For (3.16b), we evaluate first 

T 2 



T Y1 y k (A(^ n ,k)-A(iw n ,k)J =-— dkk[n B {k) + n F {k)} =-—. 



(3.20) 



Then using (3.18), we obtain 

e 4 T 2 d j f°° dp 



E R2 (iu q ,d) = - e —-^ / / ^[rn F (-p )p m 4P)A(Q-P)H^-r)n B (p )p(P)A m ^Q-P^ 



m 2 =0 

(3.21) 



where the substitutions r° — > p° and p — > q — p have been made in the second term to arrive at (3.21). 

Finally one may perform the k° and p° integrations in (3.19) and the p° integration in (3.21) by using 
the 5 functions in the spectral representations to obtain 



JkJ Prs=±1 I K P 
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and 



+ 



[A(kr + ps, k + p) A(ii/ q — kr — ps, q — k — p) + A(ps — kr, p — k)A(iv q — ps, q — p)] 

- k) • 7] A(ps 



nF{k)riF{p) Xl , . , ,s i . , , , x * /. x 

H ; [(ps - /cr) 7o - (p - k) ■ 7J A(ps - fcr, p - k) A(iv q -ps,q- p) 



kp 

n B (k)n B (p) 



kp 



(ps-fo - P • 7) + 



kp 



[{iu q - ps - /cr) 7o + (q - p - k) • 7] 



A(iu q - ps, q - p) A(ii/ g - - fcr, q - p - k) | + 



(3.22) 



S R2 (^ ? ,q) = - 



e 4 r 2 a 

4 dm 2 



I 12 \ ~ sp ^° _ (<1 _ P) ■ 7 (ii/g - sp, q - p) — ^ 
• / p s =±i I P 



-(s^ P 7o - P • 7) A (^g - «-Ep, q - P) 



np(^ 



+ ... , 



(3.23) 



m 2 =0 



where the ellipsis in each case denotes the terms with fewer than two statistical factors, which are subleading 
in the high temperature limit. We again postpone to the next section the detailed evaluation of (3.22) and 
(3.23) in Feynman gauge for a fermion at rest. 



C. The Crossed Photon Diagram 

The third and final two-loop diagram is the crossed photon diagram of Fig. 2, which topologically 
has the form of a vertex correction. The contribution of this crossed diagram to the two-loop electron 
self-energy is 

E (Q) =ie 4 [ I 7q G(Q - PhMQ - P - K) lv D» a (P)G{Q - K) lp D^(K). (3.24) 

JK JP 

In imaginary time, using the spectral representation of the propagators, this becomes 

£ c («Vq) = -e 4 f [ (n f M ) 7Q p F (P 3 ) 7 ^ F (P 5 ) 7 ,p-(P 1 )p F (P 2 ) 7 ^A(P4) i^eiit^i > 

ikj P W J Z7T J \ n,mi=l P i lUJi J 



n,m j=l 

(3.25) 



where P 1 = P, P 2 = Q - K, P 3 = Q - P, P 4 = K and P 5 = Q - P - K with again Q = (w/„q), 
P = (iOJn, p), and -ftT = (iw m , k). The Matsubara sums performed with Gaudin's method gives for the last 
factor of (3.25), 



5 1 
T 2 VTT - = 

n.m i=l r % 



nB(-Pi)n F (P2) n B (-pl)n F {p c 3 



{p\+pl-iv q ){pl+pl-iv q )[ pl+p\-pl Pl+Pl-Pl _ 
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FIG. 5: The Gaudin tree graphs corresponding to the crossed photon diagram of Fig. 2. A tree is represented by 
thick lines. As in Figs. 3 and 4 the two thin lines carry the two independent Matsubara frequencies to be summed, 
while the thick lines belong to the Gaudin tree and contain no Matsubara frequencies. 
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{Pi+Pz-iv q ){P2+Pl 



IV „ 



n F {p Q 2 )n F (pl) n B {-Pi)n B (-Pi) 
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P2 ~ P% + Wq 
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Pi + pi + pi 



+ 
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Pl+Pi+Pi 
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,0 



,0 



IV a 



nB{-Px)n F {pl) 



n B {pl)n F {pf) 



P2+P3-P5- iv q 



(P2 ~P\- Pl)i.Pl +Pl~ iv q) (P2 +P4- iv q)iPz ~ Pi ~ P5) 

nF{pf)n F {pl) n F (-pl)n F (pp 

{Pa +P5- Pl)(Pi +P3- iu q) (Pi +P5- P2)(P4 +P2~ iu q). 



(3.26) 



where the corresponding eight Gaudin tree graphs are depicted in Fig. 5. 

Substituting (3.26) into (3.25), the integrations over the three p® which do not involve the statistical 
factors or np in each of the eight terms can be performed using the definition, Eqs. (2.7a) and (2.7b), 
simply undoing the spectral decomposition in favor of the Matsubara propagator (2.4). The first and 
second of the eight terms are identical after a relabeling of momenta, while the fifth and sixth terms and 
the seventh and eight terms can also be grouped together. Suppressing the superscript on the timelike 
integration variables, we obtain 

4 f f f n f d Pl f d P2 



k J P 



■P2,k) 



+ 



+ 



• >T / . ? - n B{-pi)n F {p2) x 

j a A F (iv q -pt, q - p)7 At A F (p 2 - Pi, q - P - k)^p ua (p 1 ,p)p F (p 2 , q - k)-ypA^(ii> q 
dp2 f dp 3 

27 J ar^faWfo) x 
iaPF(p3, q - p)7m a H^<? - P2 - P3, q 

dpi f dp 4 

2tT / ^n nB i~ Pl > nB i~ P4 > x 
1 a A F (iv q -pi,q- p)7 /Lt Air(zV g - pi -p 4 ,q - P - k)j u p ua (pi, p)A F (iv q - p 4 , k^/A^, k) 
dpi f dp 5 



, v)pf{P2, q - k)7 (3 A /3Ai (^ 9 - p 2 , k) 



+ 



2vr 



2tt 



nB(-Pi)n F {-p 5 ) x 



j a A F (iu q -pi,q-p)7AiP^(P5,q- P - k)j u p ua (p 1 ,p)A F (p 1 + p 5 , q - k)jpA^(iu q - pi - p$,k) 



dpi 
2ir 



dp5 
2vr 



nB(P4)n F (p 5 ) x 
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7aA F (p 4 + Ps, q - T>htipF(P5, q - P - k)~f u A ua (iv q - p 4 - p 5 , p)A F (iv q - p 4 , q - k)7 /3 /^ t (p 4 , k) 
+ /^/^ nFfe)nF(P5)X 

J a pF(P3, q - p)7m/°HP5, q - P - k)7 !/ A" a (i^g - p 3 , p)A F (iu q + p 5 - p 3 , q - k)7 (3 A /3M (p 3 - p 5 , k) 
f dp 2 f dp 5 

7 Q A F (w ? + p 5 - pa, q - p)7^Pf(P5, q - P - k^A" ^ - P5. p)pf(P2, q - k)7 /3 A^(iz/ g - p 2 , k)} . 

(3.27) 

Expressing the fermionic and gauge Euclidean propagators A^ and A^ v , and the spectral densities p F and 
piiu j n t erms Q f ^he bosonic forms A and p defined in (2.10) and (2.9), and evaluating the expression (3.28) 
in the £ = 1 Feynman gauge, we get the covariant form for the crossed photon diagram, 

EoK,q)L =1 = -e 4 / / r ^ r d ^P{K)p{P) x 



4 JkJpJ-oo 27T J-oc 

|2n F (A; )n B (-p )C(^; Q-P;K- P)A(Q - P)A(Q - K)A{K - P) 

+ n F (k°)n F (p°)C(K; P-K + P- Q)A(Q - P)A(Q - K)A{K + P - Q) 
+ n B (-fc )n B (-p )C(Q -K-Q-P-Q-P- K)A(Q - P)A(Q - K)A(K + P - Q) 
+ [n F (fc°)n B (-p°) - n F (fc°)n B (p )] C(if + P; Q - P; K)A{K + P)A(Q - P)A(Q - K - P) 
+ [n F (k°)n F {p°) - n F (k°)n F (-p°)\ C{Q - P + K; P; K)A{Q - P)A(Q - P + K)A(K - P)} , (3.28) 

where we have defined the factor arising from the Dirac algebra, 

C(P 2 ;P 3 ;P 5 ) = C(p§,p 2 ;p§,p 3 ;pg,P5) = 7«f 3 7^f 5 7 a f 2 7 M = -8(P 2 ■ W 5 - ( 3 - 29 ) 

and relabeled p 2 — )■ /c°,pi — > p°, k — >■ — k + q in the first term of (3.27), p 2 — > k°,p 3 — > p°, k — > — k + q, p — > 
p + q in the second term, p 4 — > k°,pi — > p° in the third term, ^5 — >■ k°,p\ — >■ p°, k — > — k+q— p in the fourth 
term, p^ — > k°,p4 — > p°, k — > p, p — > — k+q— p in the fifth term, p$ — > k°,p 3 — > p°, k — > — k+p, p — > —p+q 
in the sixth term, and p§ — > k°,p2 — > p°, k — > — p + q, p — > p — k in the seventh term of (3.27), in order to 
arrive at (3.28). 

Proceeding as in the previous diagram, for each term in (3.28) the remaining two frequency integrals 
over A; and p° are performed next, using the Dirac S functions in the remaining two spectral functions. 
Then retaining only those terms which have two statistical distribution factors, we obtain 

So(*f?,q)| c=1 = ~ e4 J k J | - nF ^^ B — C(kr, k; iv q - ps, q - p; kr - ps, k - p) x 

A{iv q — ps; q — p)A(iv q — kr, q — k)A(/cr — ps, k — p) 
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nF{k)nF{p) „ , . . 

H 6 (At, k; ps, p; kr + ps - %v q , k + p - q) x 

Akp 

A(ii/ q — ps, q — p)A(zz/g — kr, q — k)A(Ar + ps — iv q , k + p — q) 

nB(k)nF(p) . , , s 

H — C (zf ? — kr, q — k; zz/g — ps, q — p; iv q — ps — kr, q — p — k) x 

A{iu q - ps, q - p) A(ii/ g - At, q - k)A(Ar + ps - iv q , k + p - q) 

nF(k)nB(p) , , , s 

— C (At + ps, k + p; zz/„ - ps, q - p; At, k) x 

2kp 

A(£t + ps, k + p)A(iz^ — ps, q — p) A(iz/ 9 — kr — ps, q — k — p) 

nF(k)riF(p) „,. , , , , x 

H ^7 C Ui/„ - ps + At, q - p + k; ps, p; At, k) x 

2kp 

A{iv q - ps, q - p) A(iu q - ps + At, q - p + k)A(Ar - ps, k - p) j + . . . , (3.30) 

where the ellipsis denotes terms with fewer than two statistical factors, which are subleading in the high 
temperature expansion. We turn next to the explicit extraction of the leading high temperature behavior 
of all the two-loop terms in the self-energy. 



IV. HIGH TEMPERATURE BEHAVIOR OF THE TWO-LOOP INTEGRALS 

In this section we analyze the integrals which result from the evaluation of the three two-loop diagrams of 
the previous section, extracting the leading order terms in the high temperature expansion T ' S> u q ,M,m, 
for fermions at rest with respect to the plasma (i.e. q = 0), and perform the analytic continuation to 
real fermion energies iv q — > u + i0 + in order to obtain the real and imaginary parts of the self-energy E. 
The behavior of the high temperature asymptotics for q = is sufficient for our purposes of analyzing 
the systematics of the two-loop self-energy. When q = only the 70 dependent part of the self-energy is 
non-zero. In this case the integrals over spatial momentum also simplify to give 

r r 1 poo poo pi 

where k = |k|, p = |p| and x = cos 6 = k • p/(kp). 



A. The Bubble Diagram 

We start with S B i(ifq) defined by (3.7)-(3.9) by evaluating the function B for q = and performing the 
p° integral with the help of the 5 function from the spectral density Pm 2 (p° >p)- The leading temperature 
dependence comes from those terms containing a product of two thermal distributions evaluated at positive 
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arguments, i.e., substituting the identities (2.20), we may discard the zero temperature 9(±E) terms and 
retain only the terms n F (\k\) and ub{\E p \) with E p = 1/p 2 + M 2 . In this way one obtains from (3.9a) 

e 4 d f°° f°° f 1 \ [ 

E B1 (ti/ 9 ,q = 0) = ^m^ — J^ dkj^ dpJ^dxkVn F (k)J2^ 



riB{E p ) 2{iv q (k + rsEp) — k(sE p — rpx) — M 2 ) — p(rp — sE p x) 



+ ..., (4.2) 



E p [2k(px - rsEp) - M 2 ] [2w q sE p - {iv q ) 2 - M 2 ] 

where the ellipsis denotes terms with fewer than two statistical factors, which are subleading in the high 
temperature limit. 

To exhibit the singular behavior of (4.2) as the regulator M 2 is removed, one may take the derivative 
with respect to M 2 inside the integral, and naively take the limit M 2 — > 0. After some straightforward 
but tedious algebra one obtains in this way, 



^4 roc roc rl f / p2 ^ \ dn B (p) 

r 7o / dkkn F (k) / dp I dx — — — ^< I h 



7T 4 Jo Jo J-i [{iv q ) 2 - 4:p 2 ] { \iu g x 2 -lj dp 



+n B {p) 



(4.3) 



X 2 — 1 \iv q P V X 2 — 1 J J (iv q ) 2 — Ap 2 \ X 2 — 1 / 

which factorizes into a hard k integral and a soft p integral. This factorization property will persist even 
when the M 2 regulator is retained, for the leading order terms in the high temperature expansion. However, 
in several of the terms in (4.3) the soft p integration is infrared divergent at p = 0, and is multiplied by 
an additional collinear divergence from the x = cos# integral at x = ±1. Thus the most severe type 
of divergence in (4.3) is the product of an infrared and a collinear divergence, each of which is linear. 
Hence one may expect generically a quadratic M~ 2 divergence from this singular behavior of the bubble 
diagram as M — > 0. This is the reason for introducing a photon mass into the calculation to regulate the 
double pole singularity of the bubble diagram. Since no infrared or collinear divergences are encountered 
in the one-loop self-energy, the divergence in the two-loop bubble diagram is the clearest indication that a 
systematic resummation of the perturbative series is needed. This HTL resummation will eventually dress 
the photon, giving it a Debye mass of order eT, and alter the counting of powers of the coupling constant. 
Retaining the mass regulator in S B i and performing the sum over r, s = ±1 in (4.2) one finds 

q4 a ( roc rl roo „2 



e 4 d ( f°° f 1 f°° p 2 

S B1 K) = ^-70^1^ dkkn F (k)j^dxj o dp-n B {Ep) 

2{iv q ) 2 Ep{px - Ep) - 2pE p (p - E p x)(px - E p ) + M 2 (iv q f + M 4 



, (4-4) 

M 2 =0 



(E 2 - vDKpx - E p ) 2 - M*/(4k 2 )} 

e appearance of kn F (k) in the firs' 
k integration is hard, i.e. dominated by k ~ T. Since M <C T, the M 4 /k' 2 ~ M 4 /T 2 <^ M 2 term in 



where ji + = [(iv q ) 2 + M 2 \/{2w q ). The appearance of kn F (k) in the first integral shows clearly that the 

r4/i,2 Aj4/rr2 ~ aj2 
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the denominator can be neglected. Then the hard k integral factorizes completely. Multiplying both the 
numerator and the denominator by (px + E p ) 2 , and observing that all terms odd in x vanish by symmetry, 
we may perform the x integral exactly and obtain 

p 2 n B (E p ) 



S B l(^g) 



e 4 r 2 a 

247r 2 i^ 70 dM 2 



[ (i "«» 2 + M2 ]/ 



+2 r dpp * Er ^4-2i V+ 1 

Jo % - Jo 



- n B (E p ) 

e 2 -n 



Ep+p 



E n 



'(I '-J": - i r : 

As we shall show explicitly in Sec. V this is exactly the form of one term which will be reproduced by a 
single insertion of a hard HTL photon self-energy in the one-loop S. 

We note that M 2 in the square bracket of the first term in (4.5) comes from the last M 4 term of the 
numerator in (4.4), which should not be neglected since it will contribute after differentiation with respect 
to M 2 . The reason for this is the collinear singularity in the x-integral which is regulated by M 2 to give 

1 2 



dx- 



(4.6) 



-i 



' {px - E p ) 2 M 2 ' 
This changes the M 4 dependence of this term in (4.5) to M 2 . 

The soft p integrals remaining in each of the three terms in (4.5) are analyzed in Appendix B. For the 
first integral we obtain from (B23), in the high temperature limit, T 3> M,T 3> oj, 



e 4 T 2 
24vr 2 w 



7o 



irT 
~M 



+ In 



M\ 
u J 



In 



-) - 2m- + ... 

CO I UJ 



(4.7) 



upon analytical continuation iu q — > oj to real frequencies. The ellipsis in (4.7) refers to terms which first 
vanish in the asymptotic high temperature limit T/M,T/uj — > oo (with M/oj fixed), and also remain finite 
in the second limit in which the regulator is removed, M/oj — > 0. In other words, terms which are divergent 
as M — > but subleading in the high temperature expansion have been discarded. It is this particular 
sequence of limits which is relevant to the comparison to the physics of hard thermal loops, and the fermion 
self-energy on-shell, in which oo ~ M are both of order of the soft scale eT, and eCl. 

The second integral of (4.5) gives after differentiation with respect to M 2 a finite contribution as M — > 0. 
Using (B24) one obtains only the purely imaginary part of the self-energy 



lj>3 



+ ... 



(4.8) 



487TW 2 

in the high temperature expansion. 

Lastly the asymptotic behavior of the third integral in (4.5) at high T is studied in Appendix B, and 
(B26) yields 



e 4 T 2 



24vr 2 



-7o 



OJ 



2M 2 



ln- + ln(4vr)- 7s j+^ 
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1 , T / T M\ / /A , \ T IT 

+-ln In -21n + In 4tt - 7i? In — - - In - + . . . 
2 u \ lo lo J \ ) M 2 uj 



(4.9) 



where is Euler's constant. 

As anticipated by the discussion following (4.3), both the real and imaginary parts of the last contribu- 
tion to S B diverge quadratically as M — > 0. This severe quadratic infrared/collinear divergence will turn 
out to cancel against other terms in S B to be evaluated below. The most interesting term is the linearly 
divergent 1/M term in (4.7), both because of its non-analyticity in M 2 and the fact that it gives rise to 
the dominant term in the real part of E, if we anticipate that after HTL resummation M will be of order 
eT. The existence of this term implies that the NLO correction to the real part of the fermion self-energy 
will finally be of order e 2 T, rather than e 3 Tln(l/e), as one might have supposed from the real part of the 
NLO one loop self-energy in (2.23). In other words, the form of the divergent terms as M — > in the high 
temperature expansion of the two-loopbubble self-energy diagram, which are not encountered at one loop, 
indicate not only the need for resummation of higher loop diagrams in a systematic way, but moreover that 
these resummed higher loop diagrams will actually give the dominant contribution to the NLO behavior 
of the real part of fermion self-energy (e 2 T), compared to the NLO one loop e 3 Tln(l/e) behavior. 

Proceeding to the second term S B2 of (3.8) one evaluates B(rk, k; rk + po, k + p; iv q — po, q — p) for 
q = with the definition in (3.7) and obtains 



g4 Q roc roo roo rl 

s B2 k) = -^70^, / o dkj o dp j_j po jj X kv j: 



n F (rk)n F (p + rk)e(p + rk)5{(p + rk) - (k + p) ) _ _ ^— _ ^ )2] J . (4.10) 

In this form one can see that because there is no k in the denominator, the k integral is clearly hard, since 
even taking account of the remaining delta function it would be UV quadratically divergent if not for the 
statistical factors. Thus k ~ T and the k integral will give a factor of T 2 . On the other hand the p,po 
integrals are soft due to the appearance of p,po in the denominators, and after differentiation with respect 
to M 2 are rendered UV finite even without making use of the statistical factors. Hence p and po will be of 
order u q or M and all expressions can be expanded in po/k,p/k. In the limit po,P *C k the arguments of 
the Dirac S are (po ± k) 2 — (k + p) 2 ~ ±2k(po =F px) and after the change of variable x — > —x in the term 
coming from r = — 1 one obtains 

e 4 d f°° f°° 1 
£ B2 K) ^ -7^4^970^72 / d Po I d PP 2 iv2 ^2 777. ^21 x 

• (4-11) 

M 2 =0 



2^ dM 2 J_ 00 J " [E 2 -p 2 ]\p 2 -^ q - Po ) 2 } 

/l f'OO r -. 

dx I dkk 2 UF(k)np(k -\- po) + (1 — rip (&))(! — np{k — po)) 5(po — px) 
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Next one can make use of a relation for the product of two Fermi-Dirac distributions in terms of products 
of a Bose-Einstein and a Fermi-Dirac distribution functions (see e.g. Eq.(3.20) of [30]): 



n F {x)n F {x + y) = (1 + n B (y))n F (x + y) - n F (x)n B (y). 



(4.12) 



Using this relation with x = k,y = po and x = k,y = —po, respectively, one obtains 

n F (k)n F (k+p ) + (1 - n F (k))(l - n F (k - p )) = -n B (-po) [n F (k + p ) - n F (k - p )] + 1. (4.13) 



Discarding the vacuum piece, expanding the two terms in the square brackets above for k 3> po and keeping 
the first non- vanishing term one substitutes (4.13) and 

, 2 dn F (k) _ f°° „ , _ ir 2 T 2 
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Jo 



dk k 



dk 



roc n i 

-2 / dkkn F (k) = 

Jo 1 



(4.14) 



into (4.11) to obtain 

S B2 (^ g ) = - 



e 4 T 2 



6vr 2 



^dM^L dX J dpXP L^^-p^-^-prf) 



M 2 =0 



. (4.15) 



As we shall show in Sec. V this form is also exactly reproduced by an HTL insertion of a hard HTL photon 
self-energy in the one- loop E. 

If one performs the po integral in (4.15), keeping only the thermal contribution which is separated with 
the help of (2.20), and taking into account the sign function coming from (2.20), one obtains 

= 4 T 2 d f 1 f 00 xp 3 n B (xp) ^ iv q 

8M 2 ./n Jo P [p 2 (x 2 -l)-M 2 } 2^ [iv q -rp(l + x)][iv q + rp(l-x)] a/ „ (j 

(4°16) 



S B2 K) = -^70 



Using partial fractioning of the denominator, changing the integration variable ptop/x and then performing 
the x integral and sum over r = ±1, one obtains 



e 4 T 2 



d | 7T 



9M 2 \ 2 Jo - - 
°° pn B (p) ( Ap 2 



p 2 n B (p) (M 2 
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pn B (p) 
p 2 - n 2 + 

p 2 n B (p) (2p 

dp^ ^-arccot 

p 2 - [i% 



, (4.17) 

io v - \"qJ ) m 2 =o 

where //± = {iiv q ) 2 ± M 2 )/(2iv q ). The asymptotic analysis of the above integrals is given in Appendix B. 
Using (B15), (B16) and the series expansion given in (B17) one obtains after analytical continuation 



S B2 (w) ~ 
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(4.18) 
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were we have again retained only the terms in the high temperature expansion which contribute at the 
order of interest or are singular in the M 2 — > limit. 

For the last term S B3 of (3.8) and (3.9c) one evaluates B(ko, k; iv q + k§ — po, k + p;po, — p) and B(ko + 
Po — iv q , k; ko, k + p;po, — p) with the definition (3.7), and obtains then 



d 



2^ l0 dM 2 

px (2rk + iv q 



7 7 7-oo 7-i - (po - iv q Y\ ^ 



PoJ 



2/cpo + r P 2 



n F (rk) n F (-p ) 



[2kr(prx + p - iv q ) + iv q {2p - iv q )] 
(r|k + p| + po — iu g )(2r\h + p|po + kpx + p 2 ) + rkpx\k + p| 



np(r|k + p| 



n F (po) 



, (4.19) 

M 2 =0 



A; 2 - (r|k + p| + p - ^<?) 2 

where |k + p| = \Jk 2 +p 2 + 2kpx. Since the large k behavior of the integrand (without the statistical 
factors) in both terms is k 2 , the k integral is hard, and k ~ T, while po and p are soft, of order v q or M. 
Expanding the expression in the curly bracket of (4.19) for po,p <C k, redefining x — > — x in the terms with 
r = — 1 and doing the sum over r yields 
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Note, that due to the relations (2.15) the first term in (4.20) is a vacuum piece and the second term 
contains only one statistical factor and can not contribute at leading order in the high temperature 
expansion. Disregarding these terms and using (4.14) one obtains 
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Performing the x and po integrals in £ #3 results in 
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, , np(p) ( 2p 
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(4.22) 

M 2 =0 



where fi± = {{iv q ) 2 ± M 2 ) / \2iv q ) . This form of the third contribution to Y*b is also obtained directly from 
the the insertion of a hard HTL photon self-energy in the one-loop £, as we shall show explicitly in Sec. V. 

The asymptotic analysis of the above integrals is given in Appendix B. Using (B15) and the series 
expansion given in (B17) one obtains after analytical continuation to real frequencies, 
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Here we have again retained only the terms in the high temperature expansion which contribute at the 
order of interest. The imaginary part neglected in (4.23) is of order e 3 Tln(l/e) and subleading to the NLO 
terms of magnitude e 2 T in which we are interested. 

The sum of the contributions (4.7), (4.8) and (4.9) from £ B1 , (4.18) from S B2 , and (4.23) from £#3 is 
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(4.24) 



Remarkably, all the M~ 2 and InM infrared/collinear divergences cancel in this leading order high tem- 
perature form, as well as 7^ and hi7r, but the M _1 divergence remains. This cancellation is in line with 
the findings of [15] which demonstrates the cancellation of all collinear and infrared divergences in the 
final result for the imaginary part of the vector-boson self-energy in the limit where the dilepton invariant 
mass is much bigger than the temperature. In our case, the presence of the photon double pole, which 
produces a more singular behavior that the fermionic one present in [15], results in the survival of the M _1 
divergence in the real part, coming from (4.7) in our final result given in (7.1a). The cancellation of the 
M~ 2 and InM terms can be understood by a simple argument, implicit already in the work of Ref. [31], 
which we discuss in Section VII. 

The most important implication of (4.24) is that because of its M _1 dependence, the photon mass 
regulator cannot be consistently taken to zero in the two-loop fermion self-energy. Only the bubble self- 
energy diagram contains this divergence. The second important feature of the above analysis is that all 
of the contributions to the two-loop bubble self-energy that are of the same magnitude as the NLO one- 
loopterms in (2.23), (or larger in the case of the M~ l term) are obtained from one hard and one soft 
momentum loop. The higher order terms which have been neglected are soft-soft. There are no hard-hard 
contributions to Eg at all. 

The remaining two-loop linear divergence oc M~ l as M — > 0, not present in the one-loop self-energy 
indicates the need for a systematic resummation of higher loop contributions and a reorganization of the 
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perturbative series. Moreover, if we anticipate that after resummation the photon will acquire a Debye 
mass of order eT, this M _1 term indicates that the correction to the real part of the fermion dispersion 
relation on-shell will be of order e 2 T, and not e 3 Tln(l/e), as might have been expected from the NLO 
one- loop correction to the real part (2.23). 



B. The Rainbow Diagram 



Taking q = in expressions (3.22) and (3.23) and taking for term S R2 the derivative with respect to 
m 2 and the limit m 2 = one obtains 
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(4.25) 
(4.26) 



{{ivq) 2 — 4p 2 ) 2 iu q {(iv q ) 2 — 4p 2 ) dp 

The first term and the very last term, involving {iv q — sp) in (4.25) can be handled in the following way. 
If k and p are hard momenta in these terms, then the sums over r, s either give integrands which vanish 
identically, or are odd under x — > — x, so vanish upon integration over x. Hence these terms have no 
hard-hard contribution and must contain at least one additional power of u in the numerator, relative to 
all terms in (2.23). Therefore they are subleading. 

Performing the sum over r and s one finds that the sum of the second and third terms of (4.25) is 

e 4 f°° f°° v f 1 dx 

The imaginary part of this contribution is subleading, cf. (BlOa), while the real part follows the hard-soft 
pattern, and does contribute at the same order as the NLO terms in (2.23). Moreover because of the x 
integral factor, it is collinearly divergent. However, taking into account the remaining terms in (4.25), i.e. 
the fourth term and remaining part of the fifth term involving kr, we observe that in both of these k must 
be treated as hard, while p is soft. After partial fractioning and using interchange of x — > —x freely, these 
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remaining terms may be written in the form [c.f. in Appendix B the derivation of (B12) from (Bll)] 
e 4 f°° „ , r . f°° , , , f 1 dx T p 2p 
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(4.28) 

In view of (B12) and (B14) the real part of this expression cancels at leading order in the high temperature 
expansion the real part of (4.27), while its imaginary part vanishes at this order. Hence combining all 
terms of (4.25) gives no collinear divergence for either the real or imaginary parts of the self-energy, and 
no contribution at the same NLO as (2.23). 

For the terms in (4.26), the factorization into a hard and soft momentum integral is again obvious. 
Performing an integration by parts in the last term of (4.26), one can rewrite (4.26) in the form 

e A T 2 iis n f°° 

dp— 

Differentiating the integrals I± (ia) defined in Appendix B gives for this contribution 
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After analytic continuation iv q — >■ uj + i0 + the R2 term becomes purely imaginary, while S R1 does not 
contribute at all at this order comparable to the NLO one-loop self-energy. Hence the contribution to the 
fermion self-energy of the rainbow diagram is simply 

EaM ~ S R2 ( W ) ~ -^^270- (4.31) 

at the same order as the NLO imaginary part of (2.23). This leading order contribution of the two- 
loop rainbow diagram arises explicitly from one hard and one soft loop integral. There are no hard-hard 
contributions to the rainbow diagram. 



C. The Crossed Photon Diagram 



Evaluating the crossed photon diagram self-energy (3.30) at q = 0, by using the definition (3.29) one 
obtains 
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In the last term we have changed x — > —x and r — > —r. Performing partial fractioning in each term one 
observes that in some cases the x-dependence partially simplifies between the numerator and the denomi- 
nator. With simple algebraic manipulations one can separate in the integrand x-independent contributions 
and contributions which are explicitly collinearly divergent, obtaining 
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(4.33) 
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Note that whenever possible we have exploited the interchanging of p ■<->■ k and x — > —x. Note, also, that 
we have catalogued all the terms as they arise in order from the previous expression (4.32), although 
actually the last term in the fourth line with two Fermi-Dirac statistical factors exactly cancels the term 
in the ninth line of (4.33). Also noteworthy is that the terms with (x 2 — 1) _1 factors, which are collinearly 
logarithmically divergent cancel whenever they multiply spatial momentum integrals which factorize into 
separate integrals over k and p, namely the last terms of the second line and the sixth line respectively. 
We shall show that the leading order terms in the high temperature expansion of (4.33) come from 
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precisely the remaining x-independent factorizable terms. Performing the interchange k -H- p in the first 
term inside the large round parentheses of the second line of (4.33) one can combine these terms to obtain 
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These leading factorizable terms contain no remaining collinear divergences. Note, that we obtained 
the product of the same integrals which already have appeared in the LO and NLO expressions of the 
HTL approximation to the 1-loop fermion self-energy (2.19). The first integral over k in (4.34) is clearly 
dominated by hard k ~ T, and given by (2.21), while the second integral is given in Appendix B and is 
soft. Through the analytical continuation iv q — > lo + i0 + , using (B6) and (B7), one obtains from these 
factorizable hard/soft terms 
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which contains both real and imaginary parts. 

It remains to show that all the remaining non-factorizable terms in (4.33) are subleading to (4.35). We 
show first that these remaining terms do not contain any contributions at all when both k and p are hard. 
The remaining non-factorizable terms are of two kinds, viz. those involving the integrand 
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multiplied by the collinearly divergent integral dx/(l — x 2 ) and either riB(p)n F (k) or n F (p)n F (k), and 
the following four terms, involving sums over r and s, 
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For the first kind of terms in (4.33), namely those involving a difference like (4.36), it is clear that there is 
no hard-hard contribution, since the integrand (4.36) vanishes for k S> to or p » uj. In fact, the integral 
over k for these terms is dominated by soft k for which we can replace the Fermi-Dirac distribution n F (k) 
by 1/2. Then the integral can be performed with the result 
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The real part of this integrand vanishes for iv q — > lo — > 0, but is proportional to 1/p for large p, the final 
integration over p with a second np(p) gives a contribution to S c proportional to e 4 w ln(T/w), which is 
subleading. In the cases where the final integration over p involves the Bose-Einstein distribution tib(p), 
there is an enhancement of T/ui and the contribution of the integral is proportional to e 4 T, which is 
still subleading. Therefore the contributions of the terms involving the difference (4.36) are subleading 
compared to (2.23) and hence NNLO. 

For the second kind of terms, that is those of (4.37), we observe that they again do not contain any 
hard-hard contributions, for when both k and p are hard, then one can ignore iv q — > oj relative to k and 
p in factors of the denominators in (4.37). In that limit, the third term vanishes upon summing over r, s, 
while the first and fourth terms of (4.37) are explicitly odd under x — > —x, and vanish under integration 
over x. By symmetrizing the remaining second term of (4.37) under x — x, we find that it becomes 
proportional to 



by antisymmetry under interchange of k and p. 

All terms in (4.37) vanish in the hard-hard limit, i.e. when oj is neglected in the denominators. The 
integral over these terms must be proportional to at least one factor of oj/T, relative to (4.35) with one of 
the momenta soft. Thus all the leading contributions from the high temperature expansion of the two- loop 
crossed photon self-energy diagram are obtained from the explicitly factorizable terms in (4.33), and given 
completely by (4.35), which exhibit the pattern of one loop integration momentum hard, and the other 
soft, just as in the bubble and rainbow diagrams considered previously. 



By explicit analysis of the two-loop self-energy we have found that the leading contributions in the high 
temperature limit which should be compared with the NLO one-loop contributions arise in all two-loop 
diagrams where one loop momentum is hard (~ T) while the other is soft (~ oj ~ eT for imaginary parts 
or over the range < k < T for real parts giving rise to logarithms). On the other hand, the basic one- loop 
HTL self-energies and vertices are derived precisely by assuming the internal loop momentum is hard, 
relative to the external momenta flowing into them which are assumed soft. Hence by replacing one of the 
loops (the hard one) of each two- loop diagram by an HTL self-energy or vertex part, and evaluating the 
resulting effective one-loop soft integrals, we must obtain exactly the same result as our direct analysis of 
the leading high temperature behavior of the two-loop diagram. This reduction of two-loop diagrams with 
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V. THE HARD-SOFT PATTERN AND EFFECTIVE HTL INSERTIONS 
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FIG. 6: Reduction of the two-loop bubble diagram to effective one-loop diagrams with HTL insertions. The framed 
part of the diagram represent the hard loop which are replaced by an effective HTL self-energy or vertex part. 

one loop hard and the other soft is illustrated graphically in Figs. 6, 8 and 10. In this section we shall prove 
that the initially rather different looking expressions obtained by the insertion of one HTL self-energy or 
vertex in the one-soft-loop diagrams of Figs. 6, 8, 10 yields precisely the same results as that obtained in 
Section IV for each diagram. This furnishes a non-trivial check by a completely different method on the 
detailed analysis of the two-loop fermion self-energy in Section IV. 

The non- vanishing three- and four-point functions in the LO HTL approximation are given in Minkowski 
space by [3, 8] 
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(5.1) 



(5.2) 
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where mr = eT/s/8 is the fermion thermal mass, mo = eT/y/3 is the Debye screening mass for photons, 
d£l = d(cos9)d(j) is the measure for angular integration and V = (1, v) with v 2 = 1. 
The HTL two-electron proper vertices obey the simple Ward identities [3] 
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while the HTL two-photon polarization function is transverse: 



(5.8) 
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A. The HTL Reduced Bubble Diagram 



Fig. 6 shows the diagrams arising when one loop is hard and another is soft in the original two-loop 
bubble diagram and the hard loop is replaced by the corresponding effective HTL two and four-point 
function. The last tadpole diagram does not contribute because there is no 4-electron effective HTL 
vertex. The contribution of the one-loop diagram with the insertion of the LO HTL photon polarization 
n (5.2) is 



Sn(Q) = ~e 2 J p 7aG(Q - P)^D^(P)U^ L (P)D^(P). 



(5.9) 



Working out the Dirac structure in Feynman gauge, we obtain two terms corresponding to the two 
structures of the HTL photon polarization in (5.2) 
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Using the relations above and setting q = 0, one can write 
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where S(P) appearing here is defined by 
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and D M 2 is defined by (2.29). In imaginary time S(iu n ,p) = —iAs(ioj n ,p), (co n is bosonic) and the 
spectral reprezentation of As(iw n , p) is 
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In imaginary time, since co n is bosonic and u q — ui n is fermionic, the Matsubara sum for the second term 
in (5.10) is given by (2.13). For the first term in (5.10) we need the Matsubara sum 



3 1 1 



- PI - lUJi p\ - pi 



n B (-Pi) 

P3+P2- iv q P%+p\~ w q 



n F {p% 



(Pi +P3- i"q){P2 +P3- iv Q 



(5.13) 



40 



Q Q-P,3 
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P,A^' P,2 



V V 

FIG. 7: The diagram corresponding to the first momentum integral of (5.10) when it is decomposed into Gaudin tree 
graphs. The dashed (solid) line means bosonic (fermionic) propagator. See Fig. 3 for other conventions on the lines. 



which was evaluated knowing that u)i,L)2 are bosonic while W3 is fermionic. The corresponding Gaudin 
tree graphs are shown in Fig. 7. 

Using in (5.10) the spectral representation of the propagators and the result of the Matsubara sums 
given in (2.13) and (5.13), one obtains 



^n{w q ) = Sni(ivg) + £112 + £113(^5) , 



(5.14) 



with 



e 4 T 2 







67r 2 " 7o dM 2 



dpp' 



g-p0 + IVqV 



00 

-00 
X 



dp e(p 



$(Po - El)n B {-po) 



(iu„ 



Po) 



£112 



e 4 T 2 
6vr 2 







dM 2 



clx 

1 PO - Xp± ) M 2 =Q 
1 /*oo ^00 

dx I dp xp 3 / dpo 

1 Jo J -00 



5(po - xp)n B (-po) 



(El-pl){p 2 - (iv q 



Po) 



(5.15a) 
, (5.15b) 
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where x = cosO and the three x-dependent terms come from the first term of (5.10), while the two terms 
independent of x come from the second term of (5.10). 

For Em the comparison between the HTL expression (5.15a) and the corresponding first contribution 
to the two-loop bubble diagram £#i is immediate, since by performing the x and po integrals in (5.15a) one 
obtains exactly (4.5). For the other two contributions, £112 an d £ri3> the expressions (5.15b) and (5.15c) 
derived by one HTL insertion on a soft loop coincide exactly with the leading order contributions (4.15) 
and (4.21) of the corresponding terms in the two-loop bubble self-energy, £#2 and £^3 respectively. The 
hard k integral in all terms of £# has converted into an HTL photon self-energy insertion, with the final 
p integral remaining soft, thus proving explicitly the reduction of the two- loop bubble to the one soft loop 
illustrated in Fig. 6. 

As a technical aside, we remark that there is a graphical correspondence between the Gaudin tree graphs 
of the two-loop bubble diagram and those of the corresponding HTL reduced diagrams, namely the first 
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FIG. 8: Reduction of the two-loop rainbow diagram to effective one-loop diagrams with HTL insertions. The framed 
part of the diagram represent the hard loop which are replaced by an effective HTL self-energy or vertex part. 

four diagrams of Fig. 3 correspond to the first diagram of Fig. 7, the fifth diagram of Fig. 3 corresponds to 
the second diagram of Fig. 7 and the last two diagrams of Fig. 3 correspond to the third diagram of Fig. 7. 
The second correspondence between Sri2 and £#2 is not obvious graphically because the fifth diagram 
of Fig. 3 gives a product of two Fermi-Dirac statistical factors coming from the independent frequencies 
corresponding to the thin lines of the fermionic bubble, whereas the second diagram of Fig. 7 contains just 
one Bose-Einstein statistical factor coming from the only independent frequency in this one-loop diagram 
represented by its one thin, dashed bosonic line. The detailed correspondence of these terms is revealed 
only by the analysis described in Sec. IV between (4.10) and (4.15), and relies upon the interesting relation 
(4.13) between a product of two Fermi-Dirac statistical factors and a mixed product of one Fermi-Dirac 
and one Bose-Einstein factor. The Fermi-Dirac factor contributes to the hard-loop HTL insertion, leaving 
behind just the soft-loop Bose-Einstein factor required to match the expression in (5.15b). 

In this way the equivalence of the high temperature asymptotic contributions of the two-loop bubble 
diagram with the corresponding HTL reduced diagram in Fig. 6 is proven. 

B. The HTL Reduced Rainbow Diagram 

Fig. 8 shows the three diagrams obtained by considering a hard and a soft loop in the original two-loop 
rainbow diagram. The last tadpole diagram does not give contribution to the fermion self-energy at the 
NLO HTL order because there is no 4-electron HTL effective vertex. The tadpole diagram containing the 
two-electron-two-photon 4-point function T^J L vanishes in Feynman gauge, since cf. (5.4) one has the 
contraction, g^ u V^V u = V^V 1 = 1 - v • v = 0. 

The contribution of the remaining effective one-loop diagram with the insertion of the LO HTL fermion 
self-energy given in (5.1) reads 

S S (Q) = -e 2 ^ 7a G(P)S HTL G(P) 7/3 ^ a (Q-P). (5.16) 
In Feynman gauge the Dirac algebra gives 

7a f f f I a = -P 2 la fi a + IP ■ V-f a fj a = 2P 2 f - 4P • VP. 
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Q Q-P,3 







P,l \J P,2 

FIG. 9: The diagram corresponding to the first momentum integral of (5.17) when it is decomposed in Gaudin tree 
graphs. The dashed (solid) line means bosonic (fermionic) propagator. See Fig. 3 for other conventions on the lines. 



Using the above expression in (5.16), the integral reduces to the sum of two terms: 

p 4 t 2 r do, r ~ p 4 t 2 (9 r 

Se(Q) = — J ^ J p yS{P)D{P)D(Q-P) + i — — 2 j p G m *(P)D(Q-P) 



(5.17) 



m 2 =0 



The propagators S(P) and G m 2(P) were defined in (5.12) and (3.13). In imaginary time S(iui n , p) = 
— iAs(iuj n , p), (w„ is fermionic) and As(iuj n ,p) has the same spectral representation as As(iw n ,p) in 
(5.12). 

The second term in (5.17) is exactly the term S R2 encountered in the calculation of the two-loop rainbow 
given by (3.15b). This can be seen immediately by applying (3.20) to the first integral in (3.15b). 

Since with the direct two-loop calculation we have seen that at leading order in the high temperature 
expansion the non- vanishing contribution comes only from £ R2 , it remains to show that the first term of 
(5.17) does not contribute at this order. In imaginary time this term is 



Ssi(ii/g,q) 



e 4 r 2 r dn 

4-7T 



n / #V (p? ' pi)p5(p2 ' p2) ^' p3)r ^n^\7' (5 - i8) 

P \ i= l J Z7T J n i=l Pi ~ iUJi 



where pi = p2 = p, P3 = q — p, u>\ = 002 = w n , and 0J3 = v q — uj n . Since oj\ and 0J2 are fermionic Matsubara 
frequencies, the result of the sum is 



1 



tY ff 1 = 

„ i=l P i ~ iUJi ^1 ~ P 2 



n F(~Pi) 



+ 



n B(P° 3 ) 

(P1+P3 ~ iu q)(P2 + P3 



(5.19) 



where the three terms correspond to the Gaudin tree graphs of Fig. 9. 

Substituting (5.19) into (5.18), performing the frequency integrals, and taking q = 0, one obtains the 
three terms 



Eei(m> 3 , q = 0) = Ssia(^g) + S S i 6 (il/ 9 ) + S S i c (w, 



1) 1 



(5.20) 



where 



Ssia(^q) 
Ssife(^g) 



e 4ji2 rl ^ roc 



np(-px) 



e A T 2 
32ir 2 iv„ 



70 / dx 



1 



1 l-x 2 Jo 



dpnp{p) 



P 



P 2 -K/2) 2t 



(5.21a) 
(5.21b) 
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<D~<D-a>- 




FIG. 10: Reduction of the two-loop crossed photon diagram to effective one-loop diagrams with HTL insertions. The 
framed part of the diagram represent the hard loop which are replaced by an effective HTL vertex part. 



Ssic(w g ) 



e A T 2 
327r 2 i^, 



-7o 



dx— 
l 1 



1 



x<1 Jo 



dp 



pn B (p) 



pn B (p) _ 2 

p 2 - {iv q /2) 2 p 2 - (iv q /(l + x)f 



(5.21c) 



For Ssia one observes that in the UV the integrand is well-behaved allowing the use of np(-px) ~ 1/2. 
Then the momentum integral gives 



1 



dp- 



1 



7T 



pxy 



+ 



4iv„ 4iv, 



' In' 



l + x 



(5.22) 



After analytical continuation one sees that the imaginary part contains an odd function in x whose integral 
vanishes. The real part is collinearly divergent, but nevertheless it is subleading because it does not contain 
the factor ln(T/w). 

Using (B6) and (B14), the real part of Ssi c cancels at leading order in the high temperature expansion 
against the real part of Ejjib- As for the imaginary part, as shown in Appendix B neither S^ib nor £ei c 
acquires a contribution at leading order in the high temperature expansion [see (BlOa) and (B14)]. This 
completes the proof that the first term in (5.17) does not contribute at leading order in the high temperature 
expansion. Hence the only contribution of the reduced HTL rainbow diagram in Fig. 8 at leading order 
in the high temperature expansion is given by the second term of (5.17), which agrees exactly with the 
corresponding rainbow diagram contribution of (4.31). 



C. The HTL Reduced Crossed Photon Diagram 

The reduction of the two-loop crossed photon diagram into diagrams containing one HTL insertion 
and one soft loop is shown in Fig. 10. The last diagram in Fig. 10 does not contribute to the fermion 
self-energy at NLO HTL order because there is no 4-electron HTL effective vertex. The contribution of 
the two remaining diagrams in Fig. 10 in which one of the original loop integral is replaced by the effective 
HTL vertex given by (5.3) are equal and give 

Sr(Q) = -lie 2 J la G(Q - P)T* TL (Q - P, Q)D™{P) . (5.23) 
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In Feynman gauge (£ = 1) one has 7 a ($ — f)V a Y = 2(Q — P) • Vy and the angular and momentum 
integrals factorize: 

p 4 T 2 frfd ? f ~ 
*N> = 4 — j^cfv J p mQ-WP). (5.24) 

In imaginary time, using the spectral representation of the propagators and performing the Matsubara 
sum coming from the momentum integral of (5.24), one obtains 

v ,. , e 4 T 2 fdn f f [°° d Po f°° dr ^ n F (r ) + n B (- Po ) 

Sr(»f„q) = ^70 / ^7^7 / / V / TT-piPo, P )p ( r o, q - P : • (5.25) 

2 7 4vr Q • V J p 7_oo 2vr J ^ 2vr Po + r - «^ 

For q = one has Q • V = iv q and one frequency integral can be performed by undoing the spectral 
representation of the propagator, while the other frequency integral can be performed by the use the Dirac 
5-function in the spectral function. Retaining only the terms with two statistical factors, one obtains 

S r(^ ? ) = 1R 2 . 70 / dp— 2 7—^h{[ n B{p) ~ n F {p)] , (5.26) 

which is exactly the result of (4.34). 

Thus, the full leading order high temperature contributions of the perturbative two-loop fermion self- 
energy coincide with those of the corresponding one-loop reduced HTL diagrams, according to Figs. 6, 8, 
and 10, just as required by the hard-soft pattern exhibited in the two- loop self-energy. 



VI. GAUGE PARAMETER DEPENDENCE AT TWO LOOPS 

Having computed the full two-loop electron self-energy in the Feynman gauge £ = 1 we consider the 
dependence of the self-energy upon the gauge parameter £ at two- loops. For £ = 1 we have shown 
that all the leading order terms in the high temperature expansion of the two-loop self-energy exhibit a 
factorized hard-soft pattern. This pattern should continue to hold for £ / 1, because in Euclidean time 
the additional factor of (1 — ^)K^K V / K 2 in the photon propagator is homogeneous of order zero under 
rescaling — > \K^. Hence it neither enhances nor suppresses the loop momentum Matsubara sum and 
integration with respect to the £ = 1 case, and we expect the same hard-soft pattern to hold. Assuming 
this to be the case, we can analyze the £ dependence of the remaining soft loop integration in the reduced 
HTL description of the last section. 

One immediate effect of the factorization into one hard, one soft loop is that all £ dependence drops 
out of the internal hard thermal loop, so that it is impossible to find any quadratic dependence on £ in 
the leading high temperature terms of the two-loop self-energy, despite the fact that the original diagrams 
contain two photon propagators. The remaining £ dependence can at most be linear. In this section we 
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analyze the linear £ dependent terms in the two-loop self-energy, assuming the hard-soft pattern and HTL 
reduction to effective one-loop diagrams. 1 We note also that owing to the transversality of the HTL photon 
polarization n^J L , the HTL reduced bubble diagram has no £ parameter dependence at all. 
For the rainbow diagram the first framed self-energy correction diagram of Fig. 8 gives 



«S s (Q) = (l-Oie 2 / L 

JK . 



1 S HTL ((5 — k) 1 



K^K V 



(6.1) 



while the second framed diagram of Fig. 8 gives the tadpole contribution 



T^ L (Q,Q,Ky- 



(6.2) 



ik ■ • 

In these expressions r^ TL (Pi,P 2 ) and T^ L (Q,Q,K) are the HTL vertices defined by (5.3) and (5.4). 

For the crossed photon diagram the reduction of one of the loops into a HTL vertex yields the £ 
dependent terms 



*Er(Q) = (1 - Z)ie 2 ^ 7„^^lf TL (Q -K,Q) + T» TL (Q, Q - K)^-^ lu 



(6.3) 



from the first two framed diagrams of Fig. 10, while the third framed diagram of Fig. 10 gives no contri- 
bution, since the corresponding LO HTL 4-fermion vertex vanishes. 

Since the HTL vertex functions obey the Ward Identities (5.5), (5.6), and (5.7) , we may contract the 
K^K U vectors into these vertices to express all the £ dependent terms in terms of the HTL self-energy 
vjHtl m fog form 



*£ 2 (Q) = 5 S r (Q) + «S E (Q) + ? S r (Q) 



(#2)2 u*#- # 



- i 



[£ HTL (Q-K)-£ HTL (Q)] 



+ [£ HTL (Q- J rT)-£ HTL (Q)] 



S HTL (Q _ R) 



0-1 



+(1 - flie 2 j K -^ {£ HTL (Q + K) - 2S HTL (Q) + S HTL (Q - K)} 
= (1 " Ik W? {"^™ " ™ « 



+ ^^-^S HTL (Q-K) ' 



<8-t 



<8-t 



(6.4) 



1 Some preliminary work on the gauge parameter dependence was reported in [16], with however a differing result for the real 
part, which we calculate fully here. 
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for the sum of these two-loop £ dependent contributions. In order to find the magnitude of this two-loop 
gauge dependent self-energy, we substitute the LO HTL condition, (ft = E HTL (Q) and add the result to 
the one-loop £ dependent self-energy (2.27), to obtain 

1 



? S on . shcll (Q) = [5 Sl (g) + «S 2 (Q)] on . shcll = (l-e)ie 2 ^^^|-^-^^^-(^ + ^) 



+< 



1 



S HTL (Q — K)- 



1 



# + £ HTL (Q + K) 



(6.5) 



on-shcll 



for the sum. 

We note that the relative sign of the first two terms in (6.5) is opposite that for one-loop £ dependent 
self-energy (2.27). Following the familiar methods these two terms can be written in the forms 

„ 1 



(K 2 ) 2 *@-fC 



(6.6a) 
(6.6b) 



for q = 0, where 

Fl EE 
F 2 EE 



d 



dM 2 
d 

8M 2 ./,: 



f°° f°° dko f°° dpo 

Jo J ~2^ PM2 ( k °>^ J ^Pop(Po,-k) 



k 2 dk 



2tt 
dko 
2tT 



n B {-kp) + n F (p ) 
ko+Po- iv q 



p M 2(k ,k) n B (k ) , 



(6.7a) 
(6.7b) 



with the limit M — > (if it exists) to be taken after integration is understood. After performing the 
integrations over all variables which are not arguments of the statistical distributions, and an integration 
by parts in the fermionic term, we find that the thermal part of F\ may be expressed in the form 



Fx 



i r 

8iv q Jo 



dE d . „o . 

; [E 2 n F (E)} 



E 2 - {iv q /2) 2 dE 1 



1 



8K) 2 Jo 



dEEn B {E) 



E + iv a 



E - iv„ 



[(E + iu q /2) 2 (E-iv q /2f 



dEn B (E) 



E + iv q 



+ 



E — iv„ 



E + iv q /2 E-iv q /2\ ' 



8iv q J M VE 2 - M 2 

where the finite M regulator must be retained only in the last integral. Likewise for F 2 , 

dEn B {E) 



-If 



(6.8) 



(6.9) 



I M VE 2 - M 2 

The high temperature asymptotic forms of each of the integrals in F\ and F 2 may be obtained by the 
methods of Appendix B, with the results 



Fi = - 



1 

4u 



7rT CO 



3iirT 
2uj 



+ ... 



(6.10a) 
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F 2 = — 
4 



vrT , T 

in — 

M M 



+ ... 



(6.10b) 



after continuation to real time, iu q — > oj + i0 + . Thus each of the first two terms in (6.5) contain both 
linear and logarithmic divergences as the photon regulator mass M — > 0. These add in (6.5) but cancel 
when combined with the opposite relative sign in (2.27), to yield the finite result (2.34) for the high T 
asymptotic form of the gauge dependent part of the one-loop self-energy. 

For the third and fourth terms in (6.5) one uses the definition of S HTL given in (5.1), performs the 
Dirac algebra, using 

<&{<& - ft) ft® - = 2(Q - K) ■ V[(Q -K) 2 @-K 2 @ + Q 2 fC] + (Q - K) 2 [Q 2 f-2Q- V<&\ (6.11) 



and the Matsubara sums, to obtain for q = 0, 
and 
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nB(-ko) + n F (po) 
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(6.12) 
(6.13) 

(6.14a) 
(6.14b) 
, (6.14c) 
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ko + ro + iv q 



(6.15) 



8M 2 J 

The Matsubara sum needed to obtain (6.14d) is given by a formula analogous to (5.19). corresponding to 
the Gaudin tree graphs of Fig. 9. In order to deal with the double poles in both the bosonic and fermionic 
propagators we have made use of the massive regulator (2.29) as before, with M the bosonic mass and m 
the fermionic one, the limits M, m — >• (if they exist) to be taken at the end. 
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The integrals in (6.14) and (6.15) are similar to those we have encountered above and in the previous 
Feynman gauge analysis. The most efficient method of evaluating them is to write them as a sum of terms 
and in each term use the delta functions to evaluate all the integrals which do not involve the Bose-Einstein 
or Fermi-Dirac statistical functions, leaving this single frequency integral for last, using relations (2.15) to 
retain only the thermal parts as integrals over positive range. In this way we obtain 
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where we have set M and m to zero wherever the integrals are finite after differentiation, retaining the 
regulator masses only when they are needed. 

The high temperature asymptotic forms of the integrals (6.16) and (6.17), evaluated by the methods of 
Appendix B are 
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(6.19) 



upon analytic continuation to real time. 

Assembling the results for the asymptotic forms (6.10), (6.18) and (6.19) in (6.5), (6.6), (6.12) and 
(6.13), with oj = irif = eT/y/8, we find finally 

„2 



? S on _ sheU (Q) = (1 - £) — 7o 
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(6.20) 



where we have everywhere neglected terms of order e 2 nif ~ e 3 T subdominant to the logarithm in the real 
part and e 2 T in the imaginary part. Although all dependence on the regulator masses M and m have 
dropped out of the sum of terms in (6.5), neither the real nor imaginary parts of the gauge dependent self- 
energy vanish, even when the on-shell LO HTL condition is used. In addition to the infrared divergence 
encountered in the two-loop bubble diagram, this residual gauge parameter dependence in a physical 
quantity such as the fermion lifetime in the plasma indicates that a resummation of the perturbative series 
in e 2 is necessary in high temperature QED. Such a resummation was performed for QED in [20, 21] within 
the HTL resummation program. After HTL resummation both the electron damping rate and the electron 
mass were calculated at NLO and shown to be independent of the gauge fixing parameter £. 

VII. SUMMARY AND DISCUSSION 



We have presented a thorough analysis of the leading order contributions to the high temperature limit 
of the two-loop electron self-energy at rest with respect to the plasma (q = 0). The three two- loop self- 
energy diagrams illustrated in Fig. 2 yield (4.24), (4.31), and (4.35) at leading order for T » 00 in Feynman 
gauge, i.e., 
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(7.1a) 
(7.1b) 
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As explained in the Introduction, the photon mass M was introduced to regulate the infrared and/or 
collinear divergences in the self-energy bubble diagram. All terms which are either independent of tem- 
perature or vanish in the high temperature limit T/uj, T/M — > oo have been omitted from (7.1). Terms 
which remain finite in the subsequent limit of removing the regulator, M/uj — > have also been omitted. 

The first important observation about the leading high temperature behavior of the two-loop self-energy 
given in (7.1) is that it is more infrared singular, both as uj — > and M — > than the corresponding one-loop 
LO and NLO self-energy, 

e 2 T 2 e 2 uj T e 2 T 

Si-ioop(w,0) = -8^7-70 + ^70 hi - + 2(1-30^70, (7 - 2) 

calculated in (2.23) and (2.24) and Refs. [10, 11] in covariant £ gauges. Since to LO, uj = mj = eT/y/8, 
the NLO one-loop self-energy gives a correction to its real part of order e 3 ln(l/e), and an imaginary part 
of order e 2 T. Using the lowest order value for uj = rrif for an electron on-shell and at rest in the plasma, 
we observe that the leading order terms in the high temperature expansion at two loops, given by (7.1b) 
and (7.1c) are comparable to the NLO terms from one-loop, again e 3 ln(l/e) in its real part, and e 2 T in 
its imaginary part for this value of uj ~ eT. The higher power of e 2 in the numerator of the two-loop 
self-energy is compensated by the more singular behavior of the denominators (7.1) when uj ~ eT. The 
1/M term in the real part of the bubble contribution (7.1a) violates the rule of finite NLO contribution to 
the real part of £ of order e 3 ln(l/e), and gives a true, uncanceled infrared divergence as the photon mass 
regulator M — > in bare perturbation theory. The singular infrared behavior of (7.1), and the uncanceled 
linear divergence as M — > in (7.1a) clearly demonstrate the breakdown of the usual perturbative loop 
expansion, even for very weak coupling e« 1. 

The gauge dependent terms of the self-energy up to two-loop order on-shell are given by (6.20), viz. 



? £ or , shcll (Q) = [^i(Q) + f S 2 (Q)] or , shcll = (1 - 70 
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(7.3) 



and are also of the same order as the NLO one- loop result (7.2). This non-cancellation of the gauge 
parameter dependence even on-shell is a second indication of the breakdown of weak coupling perturbation 
theory at high temperatures, the need for resummation and eventual restructuring of the usual (bare) 
perturbative loop expansion. 

In our direct evaluation of the two-loop self-energy in bare perturbation theory, the photon mass M 
is a regulator, introduced to define the double pole terms in (3.1), and should be taken to zero. How- 
ever, anticipating the result of HTL resummation of the NLO terms, we can estimate instead that after 
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resummation M will be replaced by the photon Debye mass, so that M ~ eT finally. With this estimate 
we observe that the 1/M term in (7.1a) gives a contribution to the real part of the self-energy on-shell 
which is of order e 2 T, parametrically larger than the e 3 ln(l/e) correction of both the NLO one-loop terms 
and the other leading order two-loop terms. Thus, the electron effective thermal mass at next-to-leading 
order in the HTL approximation is of order 0(e 2 T), not of order e 3 Tln(l/e) as one might expect from the 
NLO one- loop result, and both the real and imaginary parts of the full resummed self-energy contain order 
e 2 T corrections. In fact, the resummation within the HTL effective theory has been calculated in both its 
imaginary and real parts, and a gauge parameter independent result of order 0(e 2 T) for both has been 
obtained [20, 21]. 

The origin of this result can be understood intuitively also in the following way. From (2.19) in the LO 
HTL one-loop self-energy there is a contribution proportional to the Bose-Einstein statistical factor of the 
photon of the form 

e 2 f°° „ k 



roo u 



Because the Bose-Einstein distribution behaves like T/k at small k, this hard integral converges linearly 
as k — > 0, and as we have noted is dominated by k ~ T. However at two loops the bubble diagram has an 
additional photon propagator, and photon self-energy mass insertion proportional to H nTL /k 2 ~ e 2 T 2 /6k 2 
relative to the one-loop contribution (7.4). Thus we should expect an integral with an additional power of 
e 2 T 2 /6k 2 in the integrand, and the two- loop bubble to make a contribution to £ of 

e 4 T 2 f°°dk 1 e 4 T 3 f°° dk 



f°° dk 1 e^T 6 f 

70 Jo Te^i~n^J 



12tt 2 w ,u J k e k / T - 1 12vr 2 w J k 2 
which is now dominated by soft momenta and in fact diverges linearly at small k. When this linear IR 
divergence is cut off by the introduction of a photon mass regulator M, so that k 2 is replaced by k 2 + M 2 , 
we therefore expect a non-analytic contribution of order of 

e 4 T 3 r 00 dk e 4 T 3 



12vr 2 w 7 °y o k 2 + M 2 247twM 70 ' (7 ' 5) 
which is indeed what we find in (7.1a). Notice that this estimate does not allow for a M~ 2 or InM term, 
and is consistent with earlier approximations to the bubble diagram, c.f. [31]. There is no divergence, 
logarithmic or otherwise, from the corresponding HTL insertion of a fermion self-energy and double fermion 
propagators in the rainbow diagram. Hence there is no need for any fermion mass regulator in any two-loop 
elf-energy diagram. From (7.1a) or (7.5) the order e 2 T contribution to the real part of the self-energy in the 
HTL resummed evaluation of [21] has its origin in the appearance of a non-analytic linear 1/M dependence 
on the photon Debye mass M ~ rrif ~ eT of the two- loop (and higher) bubble self-energy diagram(s). 
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Since the HTL resummation program requires a hard-soft factorization at £ loops, in which exactly 
one of the loops is soft and the rest are hard, the hard-soft pattern we have verified explicitly at I = 2 
is a necessary condition for HTL resummation method to work, and provides a non-trivial check on its 
consistency. At £ > 2 loops the HTL resummation contains a geometric series of multiple hard photon 
and electron self-energy insertions on the basic one-loop topology of Fig. 1, but in which the bare vertex 
is corrected only to one hard loop order. Thus for £ > 2 HTL resummation when expanded in bare 
perturbative diagrams contains only a subset of all self-energy diagrams with £ — 1 hard loops and one 
soft loop. To demonstrate that the HTL resummation captures all contributions of comparable order for 
uj ~ eT, one would have to show that for £ > 2 all the contributions to £ at order e 2 T arise from 

• factorized £ loop diagrams with exactly £ — 1 hard loops imbedded in one remaining soft loop; but 

• only the subset of diagrams of this kind that can be represented as multiple insertions of either the 
HTL photon or electron self-energy, with no more than the one-loop HTL vertex insertion. 

Thus the diagrams of different topology, such as a double rainbow two-loop vertex correction on the basic 
one-loop topology should give a subleading correction to the self-energy at £ = 3 and be negligible compared 
to e 2 T for an on-shell fermion at rest. The HTL resummation program essentially assumes that this holds. 

Let us emphasize again that there are no hard-hard contributions at order e 2 T in the two-loop self-energy 
as enumeration of the generic power counting rules given in [3, 8] might lead one to expect, and moreover 
that this absence of hard-hard contributions at two-loop order is in fact a necessary condition for HTL 
resummation. The soft-soft contributions to the two-loop self-energy are non-vanishing but subdominant 
to the NLO one-loop terms in S, and will enter at higher orders in the consistent expansion in e. 

Our method has been based on isolating the leading order singularities as uj — > of the asymptotics 
of the high temperature expansion, at each loop order separately, and then comparing the magnitude of 
these contributions when the fermion is on-shell at uj ~ eT/\/8. This suggests that collecting all the 
contributions which are subdominant to the NLO ones we have been considering, i.e. the NNLO order 
contributions, would require two soft loops and 1—2 strictly hard ones of order T, which would require either 
at least a three-loop computation, or an inductive general argument to verify explicitly. If this is the case, 
then collecting infrared singularities order by order in the high temperature expansion at each £ becomes 
possible, and a reordering of the usual loop expansion can be carried out in a completely systematic and 
gauge invariant way. At NLO this reordering is just one-loop HTL resummation, while at higher orders it 
may require additional resummation of both self-energy and vertex parts. 
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The HTL resummation procedure works when all energies and momenta in the remaining loops are soft 
~ eT. This one-soft-loop HTL resummed contribution should always give the order oj/T ~ e corrections 
to the leading order HTL self-energy. On the other hand, uncritical use of HTL resummation has no 
justification a priori when extended to higher orders in oj/T in the high temperature expansion, or in loop 
integrals containing hard momenta, and may lead to incomplete or incorrect results when misapplied in 
this way. This seems to be the case in the use of HTL resummed propagators in the photon polarization 
function for the calculation of the photon damping rate in QED, or photon and dilepton production in 
QCD [32-34]. In the photon self-energy (unlike the fermion self-energy considered in this paper) the NLO 
result vanishes, because of the appearance of the statistical factor 1 — 2np(uj/2) ~ uj/AT which produces a 
further suppression [33] . Because of this suppression the first non-zero contribution to the photon damping 
rate in QED and the complete calculation of the photon and dilepton production rate in QCD both require 
additional diagrams and depend upon additional physical processes (bremsstrahlung) than those captured 
in the one-loop HTL resummed photon polarization tensor *H^ U . Based on the results of this paper, we 
would conjecture that a completely analogous calculation of the two-loop perturbative photon polarization 
in QED would give a vanishing contribution from both hard- hard and hard-soft loops, and yield a non-zero 
contribution only when both loops are soft, which is down by (u/T) 2 ~ e 2 and is NNLO compared to the 
LO one-loopresult e 2 T 2 . Hence it cannot be captured by the usual HTL resummation. 

Once HTL resummation has been carried out, one obtains a new HTL improved self-energy which 
contains the correct physics to NLO, at the soft scale eT. Thus we would expect that the main results of 
our detailed two-loop study for fermions at rest q = should continue to hold for fermions which are slow 
moving, |q| < eT. However, there is no necessity for the expansion in q 2 /u 2 to be analytic at q 2 = 0. Non- 
analytic behavior has been reported in scalar QED [35, 36], but the spinor QED or QCD cases have not 
been fully investigated to our knowledge. It is known that even the HTL resummed fermion self-energy is 
inadequate to describe the damping rate or energy loss of a fast moving fermion with |q| » eT [13, 37], and 
a further or different sort of resummation is required [29] . The appearance of additional logarithmic terms 
of the form ln(q 2 /w 2 ) for |q| » eT does not mean that HTL resummation is inconsistent, but only that it 
must be properly applied. This example of the damping rate of a fast moving fermion, as those mentioned 
above of the photon damping rate in QED and both the photon production and dilepton production rates 
in a hot QCD plasma simply lie outside of the range of strict applicability of HTL resummation, which is 
valid for soft external four-momenta only. 

To investigate the breakdown of HTL resummation, the systematic method which we have used in 
this paper for the two-loop perturbative self-energy could be applied to the HTL one-loop resummed self- 
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energy (photon or fermion) as well. Thus, if this HTL resummed self-energy is given off-shell for all lo 
and q, one could examine its high temperature limit, or equivalently, its infrared singularities at q = 
and uj — > 0, as well as its behavior at large |q| S> eT. Then, depending on the result one could examine 
terms at two and higher loops in the effective HTL theory to determine when they become of the same 
magnitude as the LO HTL results. It is known that the breakdown of the HTL resummation program 
is associated with the appearance of collinear divergences [38] and requires an additional resummation 
of ladder-like diagrams corresponding to multi-loop vertex corrections not contained in the usual HTL 
resummed self-energy [39]. Such ladder sums are also required to compute transport coefficients in hot 
gauge theories [40, 41]. Quantitative control of the breakdown of the HTL resummation program by the 
extension of the methods utilized in this paper to determine the most infrared sensitive contributions at 
each HTL resummed loop order can provide detailed guidance to its proper generalization to more complete 
resummation algorithms in hot gauge theories. 

If one finds that higher order terms beyond NLO can be resummed in a consistent gauge invariant way, 
this would provide the theoretical justification for a systematic extension of the HTL program to higher 
orders, including also resummed vertex parts and ladder diagrams, which might enable the extension of 
resummation methods to lower energy and momentum scales and the hydrodynamic limit, where weak 
coupling methods are usually thought to fail entirely, particularly in QCD. 

Physically, taking account of infrared and collinear singularities by a reordering of the loop expansion 
and resummation of self-energies and proper vertices relate the dressed quasiparticles and many-body 
interactions in the medium to the single particle excitations of the theory in the absence of a medium, 
i.e. in the vacuum at zero temperature and density. The dispersion relations of the dressed quasiparticles 
for low frequencies and soft momenta are required to describe the collective modes of the plasma in the 
hydrodynamic limit, and to determine transport coefficients such as electrical conductivity, and shear and 
bulk viscosity. A thorough understanding of the quasiparticle degrees of freedom at high temperatures 
or densities and the interactions they undergo in thermodynamic equilibrium or under slight departures 
therefrom is thus an entry point into the systematic study of non-equilibrium processes in the plasma and 
self-consistent methods for studying them systematically as well. 
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Appendix A: Gaudin's Method for Matsubara sums 

The Gaudin method [22] relies on a systematic decomposition of a product of N fractions, each linearly 
depending upon Matsubara frequencies and corresponding to a finite temperature propagator, into sums 
of products of fractions, in each of which the only Matsubara frequencies which appear are the independent 
ones to be summed. Since the number of such independent Matsubara frequencies in a Feynman diagram 
with £ loops is I, no matter how many propagators (N) it may contain, the independent Matsubara sums 
are more easily performed in this method when I < N. 

The necessary systematic decomposition of the original product of denominators containing the Mat- 
subara frequencies in a given Feynman diagram can itself be represented graphically. For simplicity and 
definiteness we consider only the case of one external Euclidean frequency v as in the case of the multi-loop 
Feynman self-energy diagrams for £ treated in the text. 

To each Feynman multi-loop diagram one associates a set of Gaudin graphs Q. The set is composed of 
tree diagrams specified by all the possible ways of connecting every vertex of the original Feynman diagram 
without forming any closed loop. In each of these Gaudin graphs the propagator lines which belong to the 
tree can be drawn with thick lines and form a set denoted by Tg. The complement of this set, namely the 
propagator lines of the original Feynman diagram not contained in the given Gaudin tree graph is denoted 
by Tg- If the total number of propagator lines in the original graph is A, we must have 

N Tg +N % =N 1 (Al) 

where Nj- g is the number of lines in the Gaudin tree graph g G Q , and Nj- the number of lines in its 
complement. For an original Feynman diagram with I loops, in fact Nf- g = t. The propagator lines in T g will 
be treated as the lines with the independent Matsubara frequencies to be summed oj^, k = 1, . . . , Nj- , while 
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the frequencies in T g , denoted flj, j = 1, . . . , N-j- g are treated as dependent. The frequency Qj{iv,pi, . . . ,p^) 
of each line of T g is expressed in terms of the independent real frequencies p° k , k = 1,...,£ by simply 
replacing each iuj], with p® and the external imaginary frequency iu, using energy conservation at each 
vertex. 

The general Gaudin decomposition formula then reads 

N 



TT c = y 

1 1 o° - iuJi ^ 



e^TT 1 TT ^ k 



(A2) 



The time arguments in the exponentials are introduced as regulators so that the sums over the independent 
Matsubara frequencies in Eqs. (A3) and (A4) below converge. These time arguments are somewhat 
arbitrary and will be taken to zero at the end. They are needed only to guarantee that all sums converge, 
and results of Matsubara sums which are absolutely convergent will not depend on how the Tj are assigned. 
For definiteness one may choose the following rules to fix the r\ 

• Specify the orientation of each line of a diagram, label it with integer number (k) and associate to 
each line a positive number t\~. A possible choice is = ke, with e positive. 

• Tk for a given line is calculated by choosing a loop containing that line and adding the Tj of all the 
lines of the loop, formed with the line k £ T g and lines of the tree 7~ g , with ± sign if the orientation 
of the line agrees with/is opposite to the orientation of the line k G T g . 

• T e is determined by the r's of the lines of T g which connects the two external lines. In this case the 
reference orientation is the orientation of the external line. 

To perform the sums over the independent Matsubara frequencies set = ke, let e — > + and apply 
the summation formulas 

T E o_ • =e{T)n B {e{T)p°)e? T , for co n = 27mT, (A3) 

o_ • = -e(T)n F (e(T)p°)eP» T , for u n = (2n + l)vrT , (A4) 

valid for — /3 < r < fi. This associates a statistical factor to each line of T g . The e tuJnT factors are necessary 
as regulators so that the Matsubara sums are well defined and convergent. After the sums have been 
performed, the limit e — > 0, r — > may be taken. 

The simplification of the Gaudin method over other methods for performing the sums results from the 
fact that the frequencies of the last product in (A2) are disentangled and independent of each other. For 
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each of them the Matsubara sum may be performed separately, using only the fundamental formulas (A3) 
and ( A4) , each independent sum yielding its own statistical factor depending on whether that line is bosonic 
or fermionic. With all Matsubara sums performed, the finite temperature amplitude is expressed in terms 
of N timelike momentum integrals over p®, i = 1 , . . . , Nj- g and p®, j = 1, . . . , Nj- g and £ spatial momentum 
integrals. The first set of N — £ timelike integrations can be performed immediately by undoing the spectral 
representation of the propagators. For subtleties and a more extensive discussion of the Gaudin method 
the reader is referred to Ref. [23]. 



Appendix B: Asymptotic analysis of integrals 

We present here the high temperature expansion of various integrals which were encountered in the 
text. We start with the evaluation of the integral 

f°° 1 
/±(io) = ^ dxg±(x) rr2 _ ( ^ 2 



1 f°° 
2iaJ dX9±{x) 



x 2 — (ia) s 
1 



(Bl) 



x — la x + xa 

where g±(x) = x/(exp(x) ± 1). This integral appears in the calculation of the one-loop self-energy at NLO 
in the HTL approximation, and is typical of the various integrals over soft momenta that are encountered 
in the two-loop self-energy. Its real part emerging after analytical continuation to real frequencies 

»=i^ +i0+ < B2 > 

was evaluated only numerically in [10]. 

We note first that although w < T, we cannot set a = in the integrand, since then g_(x)/x 2 ~ x~ 2 
as x —7- and the integral would be linearly divergent. As we will see, the linear divergence survives in the 
imaginary part of I- after analytic continuation but not in the real part, giving a logarithmic dependence 
on T/lo instead. Indeed integrating (Bl) by parts twice, we obtain 

I±(ia) = ^ <7±(0) [ln(m) - ln(-m)] - ^g' ± (0) [\n(ia) + ln(-m)] 
1 f°° 

— / dxg±(x) [(x + ia) ln(x + ia) — (x — ia) ln(x — ia)\ . (B3) 

Ha Jo 

Under the analytic continuation (B2), 



ln(m) -»• In [-) , In(-io) ->• In [-) - m , (B4) 
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so that 

/ (,) . \ iir r T I 7\ tit 

+ ... , (B5) 



^(^ + i0+ ) = ^^°) + ^°)[ 1 "© + 



T \ in 
~2 



where the ellipsis denotes (real) constant terms, and terms that vanish as oj/T — > 0. 

Using then <7+(0) = 0, g~{0) = 1, and g±(0) = ±1/2, we find that to leading order in the high 
temperature expansion 

7o V-(w + i0+) 2 /4 w v y 

Note that for T » w only the integral with the Bose-Einstein statistical factor has an imaginary part. 
These results to leading logarithmic order may also be obtained more rapidly by simply substituting the 
first two terms of the Taylor expansion of g±(x) into (Bl), and cutting off the logarithmic UV divergence 
by a cutoff of order of T. 

An analytic calculation which captures the large temperature behavior of (Bl) beyond the leading order 
given in (B6) and (B7) can be performed with the standard technique presented in detail in Appendix C 
of [42]. Using the identities 

1 1 %. — ■ 



\^ E .2 , 772 7-y ( B8 ) 



exp(x) ± 1 2 ^ x 2 + oj+(nY 

v ' n=— oo =t v ' 

where w+(n) = (2n + l)ir and w_(n) = 2n7r, one separates the integral given in the first line of (Bl) 
into two integrals which are regularized by inserting x e in the integrands. Interchanging the order of the 
summation and integration one uses partial fractioning and exploits the properties of the sums to write 

r i r°° , x 1+e , r°° x e 

I±(ia) = lim < ± - / ax^ ? + t± / 



e-s>o [ 2 y o x 2 + a 2 y o x2 + fl2 

±2G Vo ^x 2 ^ 2 5 u4(n)-a 2 ^ 2 5 ojl(n)-a 2 J dX * 2 + c4(n)/' (B9) 

where Z± = — (— l±l)/2. The second term is the contribution of the static Matsubara mode n = in 
7_(a). In the second and third terms the limit e — > can be taken. In the last term one does the integral 
and then one expands the summand around a = in order to be able to do the sum over n. Note, that for 
bosons this expansion can be safely done once we separated the static mode. The limit e — > is taken at 
the very end which results in the cancellation of the 1/e type singularity present in the first and last terms 
of (B9). With these sequence of steps one obtains 

1 A„ 2a , ^\ , 7ra 7 C( 3 ) 2 , rni 3x 



J+(ia) = "2 ( ln V +7£ j + T" 8^ + ° ( )} (B1 ° a) 
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I- ha) = h - In — 

y ' 2a 2 V 2tt 



(BlOb) 



where a = v q j(2T) and 7^ is the first Stieltjes constant (Euler's constant). 

We can use the result (B5) or (BIO) for the integral (Bl) to evaluate the integral appearing in (4.28) 
and (5.21c) at leading order in the high temperature expansion. One defines the integral 
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where f-(x) = n#(x) and /+(x) = np(x). In the limit k 3> p,v q one approximates —2kp{x — rs) — 2iv q r{k+ 
rskp — rivg/2) ~ — 2k[p(x — rs) + riv q ], then one uses partial fractioning and the change x — > —x in some 
terms to rewrite the integral in the following form: 
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Then one uses (B6) and (B7) to write 
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where l± = — ( — 1 ± l)/2. The imaginary part vanishes since the integrand is odd in x. Keeping in the real 
part the terms which are leading order in the temperature one has 
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With the same method one can perform the asymptotic analysis of the integrals appearing in (4.17) and 
(4.22). Introducing a new variable of integration x = p/T, the notation a = u q /(2T),b± = (v 2 ±M 2 )/(2v q T) 
and taking the derivative with respect to M 2 of the various terms in (4.17) and (4.22) one obtains 
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(B15b) 
(B15c) 
(B15d) 
(B15e) 
(B15f) 
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where again g±(x) = x/(exp(x) ± 1). The only term where one cannot take the limit M 2 = is 
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Here we expanded in powers of M. In the expressions above I±(ia) is given by (Bl). 
Then the method described between (B8) and (BlOa) gives the following series: 
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where the constants are given in terms of the first two Stieltjes constants and 71 and the Glaisher's 
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In order to derive (4.7), (4.8), and (4.9) first we define two integrals, 
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where /i± = {{iv q ) 2 ± M 2 )/(2iv q ), E 2 = p 2 + M 2 , and ln(c) = 7^ — 1/2 - ln(47r). The integrals were 
evaluated in the high temperature expansion with the method described between (B8) and (BlOa). The 
first term in the expression of J_ (/x_) is the contribution of the static Matsubara mode and can be obtained 
also with the approximation nB(E p ) ~T/E p . 
Then using that 
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the first two integrals in (4.5) can be written as 
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Using the expressions given in (B18) and (B19) one obtains 



d 



dM 2 



(( 



lV„f +M 



roc 

Jo 
) r 

^Jo 



dp 



p 2 n B {E p ) 



E v E 2 -^\ 



dp- 







" 1 


2i ' 




~ vrT 










M ~ 


iv q 


) 


1 


ttT 






8 


l-r. — 
AlV q 





, M 1, T „ 

+ In - In — + C, 

IV a I IV„ 
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where C = (1 + je — ln(47r))/2 and the expressions on the right hand sides were expanded for small M. 
After analytical continuation (iv q — > oj + i0 + ) one obtains (4.7) and (4.8). 

For the last term of (4.5) one may introduce x = E p /T as a new variable of integration and obtain 
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where \i = M/T and b = (v 2 - M 2 )/{2v q T). In order to obtain the contribution of (B25) for M, \iv q \ <C T 
at leading order in the high-T expansion one splits the integral into two pieces by introducing the point x 
satisfying 0</i,!)<i<l. For x > x one expands the logarithm for large x and use the method described 
between (B8) and (BlOa). For x < x one expands x/{e x — 1) ~ 1 — x/2 and one scales the variables with 
\i. Writing the logarithm as a difference of two logarithms one introduces in the resulting two integrals the 
argument of the respective logarithm as a new variable of integration. The integrals can be done with the 
symbolic program Mathematica. Adding the two parts the x dependence cancels and the result is 
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where Cb = "^ — 7i + i2+5 ln 2 (47r) — (7^ + \) ln(47r). After analytical continuation one obtains (4.9). 
The high temperature asymptotics of integrals of the form 

dE 
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(B27) 



10 VE 2 - M 2 

with f(E) a regular differentiable function of E at E = 0, encountered in Sec. VI can be handled by 
similar methods. Introducing an intermediate scale E, such that M <C E «T, one can write 

' M dE_ 
i M VE 2 - M 2 

where the ellipsis gives terms of order M/T or smaller. Performing the first integral exactly and integrating 
the second by parts twice gives 
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+ — [En B (E)f(E)]ln- _ + 



irT T 
= 7777/(0) + Tf'iO) In— - 
2M J w J y ' M 




(B29) 



in the high temperature limit. This gives the results (6.10) and is used in obtaining (6.18) and (6.19). A 
result analogous to (B29) is easily obtained for the Fermi-Dirac distribution rip. Finally we note that the 
same method of introducing an intermediate scale, \w^\ <C E <C T may be used to extract the asymptotic 
behavior of the finite integrals in Eqs. (6.16) as well. 
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